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FOREWORD 


The overall objective of this NASA program has been to develop and 
implement several computer programs suitable for the design of lobe 
forced mixer nozzles. The approach consisted of extending and 
existing analytical nacelle analysis to handle two stream flows where 
one of the streams is at a higher energy. Initially the calculation 
was set up to handle a round, free mixer including satisfying the 
Kutta condition at the trailing edge of the mixer. Once developed and 
calibrated, the same analysis was extended to handle periodic 
boundary conditions associated with typical engine forced mixers. The 
extended analysis was applied to several mixer lobe shapes to predict 
the downstream vorticity generated by different lobe shapes. Data was 
taken in a simplified planar mixer model tunnel to calibrate and 
evaluate the analysis. Any discrepancies between measured secondary 
flows emanating downstream of the lobes and predicted vorticity by 
the analysis is fully reviewed and explained. The lobe analysis are 
combined with an existing 3D viscous calculation to help assess and 
explain measured lobed data. 

The program also investigated technology required to design forced 
mixer geometries for augmentor engines that can provide both the 
stealth and performance requirements of future strategic aircraft. 

For this purpose, UTC's available mixer background was used to design 
several preliminary mixer concepts for application in a exhaust 
system. Based on preliminary performance estimates using available 
correlations, two mixer configurations will be selected for further 
testing and analysis. 

The results of the program are summarized in three volumes, all under 
the global title, "Turbofan Forced Mixer Lobe Flow Modeling". The 
first volume is entitled "Part I - Experimental and Analytical 
Assessment" summarizes the basic analysis and experiment results as 
well as focuses on the physics of the lobe flow field construed form 
each phase. The second volume is entitled "Part II - Three 
Dimensional Inviscid Mixer Analysis (FLOMIX)". The third and last 
volume is entitled "Part III - Application to Augmentor Engines". 
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SECTION 1.0 
INTRODUCTION 

A forced mixer is a device that is used on gas turbine engines to internally 
mix the hot turbine efflux with the cooler, lower velocity fan bypass or 
secondary stream. The principle motivator for doing this on commercial gas 
turbine engines is to reduce the jet noise associated with the high energy 
core stream. Also, when designed properly, the mixer can also achieve an 
increase in gross thrust while realizing noise reductions. Figure 1 shows a 
typical separate and mixed flow nacelle. Note that to achieve the mixing, the 
duct surrounding the engine must be lengthened. This adds weight, as does the 
mixer and the larger centerbody plug. Normally, and particularly for short 
range aircraft applications, these weight penalties offset the thrust 
improvements. In addition, the increased difficulty of intergrating the longer 
nacelle into the aircraft flow field without incurring interference drag 
penalties, in the past has prevented launching new engines into service as 
mixed flow nacelles purely on a performance basis. However, with increasingly 
stringent noise regulations, mixers are being considered for future commercial 
engine applications. Mixers are also being considered in military applications 
as variable area devices for varying cycle match and as a way of spreading the 
hot turbine exhaust ahead of afterburning flameholders. This report presents 
an analytical method which has been developed for the mixer lobe flowfield. 
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Figure 1 Typical Nacelle Exhaust Configurations 
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Conventional or commercial forced mixer geometries consist of periodic lobe 
structure that can be described in terms of a number of geometric features. 
Figure 2 illustrates two cross-sectional views of the lobe. The scarf angle or 
lobe cutback angle is used to reduce lobe length with increasing penetration 
into the core. Scalloping is a cutout of the lobe lateral surface used to 
minimize lobe wall structural problems while promoting tangential mixing. 




Figure 2 Mixer Lobe Geometry Definition 
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For many years the mixer was designed using a trial and error experimental 
approach, wherein limited traverse and performance data was used to refine 
desiqn concepts. 1 “ 3 More recently, "benchmark” experiments *» = resorting 
to high response and LDV instrumentation have probed the mixing chamber in an 
attempt to explain the mixing process and its driving mechanisms. These 
experiments confirmed that the mixing process is a viscous dominated process, 

6 and that the primary driving mechanism is the secondary flows generated in 
the lobed region of the flow. Several researchers have proposed a variety of 
inviscid and viscid processes for producing the secondary flow, but asye , no 
attempt has been made to analytically model them. ( Anderson and Povinelli > 
have lumped these terms together under a nerenc vorticity label, 
analytically simulating its effect in terms cf a vortex sheet distributed 
along the lobe exit surface. Such an approach has been used to generate inlet 
conditions for a viscous marching analysis in the mixing duct. The results of 
these calculations, as seen in Figure 3, have been shown to rea . i s r i cal ly 
simulate observed flow mixing patterns. The purpose of this paper 1 * t0 
develop an inviscid analysis which, in conjunction with a lobe boundary layer 
analysis, can predict the flow over the mixer lobes, thereby obtaining the 
conditions needed to initiate a marching viscous calculation in the downstream 

duct. 

The analysis presented in this report is an attempt to examine one P a ^ular 
inviscid secondary vorticity generator called the flap .jorticity schenv e 
(Figure 4). In this model, the vorticity is associated with the periodic lift 
distribution produced by the periodic lobe trailing edge. The feasibility o 
an inviscid method predicting an observed level of flow 
dependent on the axial flow component remaining attached to the surface. 
Surface pressure and shed vorticity distributions will therefore be examined 
to calibrate the method. 
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Comparison of computed temperature profiles with experimental 
data. Case 12C Mixer, Tf/T c = 0.74, Station 21. 
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SECONDARY FLOW GENERATION, PASSAGE VORTEX 


Figure 4 Secondary Flow Generation Models 
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A. Overview 


SECTION 2.0 
ANALYTICAL APPROACH 


The forced mixer consists of a convoluted lobe section and a mixing chamber. 
Observations have proven that the lobe region is responsible for the secondary 
flow generation while the downstream duct region to the nozzle exit plane 
produces the resultant flow mixing. It is reasonable therefore to propose 
zonal analytical approach, wherein regions are treated using locally 
applicable techniques. The convoluted lobes can be viewed as a ring wing with 
a periodic "spanwise" loading distribution, r(e) , as is shown in Figure 5. The 
periodic lobe cambering produces a nonuniform loading distribution and a 
corresponding shed vorticity field where strength varies periodically in e. 

The shed vorticity is associated with crossflow velocity field (secondary 
flow) that "mixes" the flow as it is convectec! downstream. The vorticity is 
then stretched and eventually dissipated in the mixing chamber through the 
action of viscosity. 



Figure 5 - Ring Wing Analogue of Forced Mixer 

A complete inviscid treatment of the three-dimensional lobe region is still a 
difficult problem, due to geometrical complexity, multiple energy streams and 
compressibility. Considering the lack of certainty about the relevant driving 
mechanism for the secondary flow generation, it is appropriate to consider a 
more approximate analysis that can verify the magnitude of the flap vorticity 
model. By neglecting the effect of wall boundary layer development, the flow 
can be considered as irrotational regions separated by a vortex sheet. 
Downstream of the lobe structure, the flows can initially considered to be 
inviscid also, with the wake modeling all of the energy jump. An analysis in 
terms of a velocity potential, defined within each region, is therefore 

possible. 
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The m'xer lobe surface is characterized by two length scale ratios, which 
provide relative measures for the lobe height, axial and azimuthal variations. 


e 



MaR 
z 2 ~ THK 


where N is the number of lobes, L the axial length and a R the lobe height 
above some mean reference radius. For many current designs, the lobes are 
axially slender (ej«l) and the local Mach numbers are low enough that a 
small disturbance model may be introduced as a means for treating the salient- 
features of the lobe mixer problem. In the next section a unique small 
disturbance formulation will be developed that will analytically uncouple the 
e variation and reduce the problem to the solution of a sequence of 
axisymmetric problems. These problems can be solved by taking advantage of 
previous experience a whereby the effects of power addition and "exact" 
surface boundary conditions were modeled. In contrast to this earlier work, a 
finite volume cylindrical grid formulation similar to that of Wedan and 
South s y is used. This approach yields a straightforward treatment of 
extremely general geometries. Although Wedan's calculations have been applied 
to essentially symmetric geometries, the method presented below will have no 
such restrictions. 

The power contribution and lobe loading result in a potential -lump [rf] across 
the wake shed from the lobe trailing edge. The corresponding induced secondary 
flow has a net circulation r . By appropriately choosing the closed path of 
integration over half a lobe, as shown on Figure 6, the circulation inteqral 
reduces to 


r= 




( 1 ) 



Figure 6 - Schematic of Lobe Trailing Edge Integration Path 


6 


In the following sections, a potential analysis for the forced mixer will be 
developed. The complex three-dimensional problem previously described will be 
linearized permitting an uncoupling of its 0 dependence in the governing 
equations. Special treatment of surface and Kutta conditions however will be 
needed to insure compatibility with the governing equations and the physics of 
the problem, while also uncoupling their e dependence. Finally the treatment 
of power addition within a potential format will be developed. 

B. Potential Flow Analysis 

The inviscid lobe analysis will be applied to the flow domain between the fan 
and core flow discharge plane and a downstream plane in the mixing chamber, 
schematically shown on Figure 6. The plane should be displaced sufficiently 
from the nozzle exit plane so as to avoid non! near compressibility effects. 
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Figure 7 Mixed Flow Macelle Analysis Domain 

The governing equation of mass conservation, applied to an arbitrary control 
volume in space, yields 


<0p V . 7f dA = 0 ^2) 

The nondimensionalized mass flux is defined in terms of a perturbation flow 
fron an upstream subsonic x-"aligned" flow as follows 


p v = 


+ R T + i> i + — i i 

p x x r r r v q 0 


(3) 
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where & 2 = i - n 2 oo and i is the perturbation potential. Only the linear 
terms have been retained in Equation (3). The expansion is made assuming that 
all surfaces can be considered to have slender shapes in the streanwise 
direction. 

Separation of Variables and Numerical Approach 

Appropriate surface boundary conditions can be derived from a flow tanqency 
condition given by: 


v ' VF = 0 (4 

where 7 , the velocity vector is given in terms of a velocity potential and F 
is a general surface described in terms of the r, rf, x cylindrical coordinate 
system 


F ( r, x, e) = r -f(x,e) = 0 ( 5 ) 

If the velocity potential is assumed to be a perturbation about the upstream 
axial flow and if the perturbation velocity components are assumed small 
relative to the upstream flow, then the axial perturbation potential 
contribution can be neglected and the surface boundary condition reduces to an 
expression for the surface radial velocity on a mean surface R n (x). 

f 

0 r = f x + <t>g (6) 

m 

The perturbation or small disturbance approximation is equivalent to limiting 
surface slopes to order ej. This linearization of the boundarv condition is 
needed to render the overall problem separable, as will be seen shortly. At 
first glance, it appears that the last term in Eq. (6) is of order e i, 
however, no such perturbation restrictions has been imposed in the azimuthal 
direction. The terms, f x , f e are known functions which describe the lobe 
geometry . 

Problem closure is obtained by imposing a quasi one-dimensional analysis for 
definition of the inlet flux and a Kutta condition to uniquely set the net 
circulation. In order to simplify the analysis, a cylindrical coordinate 
system orientation was used to simplify evaluation of the flux integral (?). 
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Normally, a body conforming or sheared Cartesian grid would appear the logical 
choice for a coordinate mesh. By treating the mixer with a pure Cartesian 
grid, the following analysis will be more tractable. Furthermore, one can also 
substitute storage of large arrays with the additional complexity of irregular 
boundary /mesh intersections. Analysis of complex geometries in two- and 
three-dimensions can be easily treated using such a scheme. The mixer lobe 
geometry is assumed to have no scalloping and scarfing so that the trailing 
edge will align with the mesh. Extensions of the analysis to include scarf 
effects will be discussed in an Appendix. Consider the situation where the 
general three-dimensional lobed mixer problem is perturbed about some mean 
surface. It is possible to avoid analyzing the full three-dimensional problem 
be recognizing that the flow is a periodic function of the number of lobes. 
Solutions to Equation (2) and any appropriate boundary conditions are 
therefore assumed to be separable, i. e., 


i ( x , r , e) = g(x,r)h(e) 


(7) 


Combining Equations (2,3,7), the governing flux balance equation terms can be 
appropriately separated into terms that are either a function of (r,x) or of e 
alone. 


E-W 





■g rdr + J rg dr 
* H-$ r 



dr d x 
r 



( 8 ) 


where K is the separation constant and the E-W, N-S and aA integrals are 
elemental areas evaluated in the (x,r) plane. The e component of Equation (8) 
can be recast to identify the periodic nature of the separated variable. The 
solution for h(e) includes linear combinations of trigonometric functions, 
where appropriate application of a symmetry boundary condition at the lobe 
crest simplifies this to 


h 


k " 



k = 0,1,2, ...MH 


(9) 


where K = ke/e 0 , e 0 is the half angle of the lobe (crest to trough), MH is 
the number of harmonics used in the Fourier series, and is a sequence of 
unknown coefficients still to be determined. The k = 0 solution corresponds to 
the axi symmetric solution limit. Since the separation constant can take on 
multiple values, the potential assumes a more general form 

6 - Z 9 k \ 

k 
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The axi symmetric component of Equation (8) can be evaluated for an arbitrary 
point in the flow field on integrals aligned to the cylindrical mesh (Figure 
8a) to produce 




1 g. drdx = 0 
r k 


( 10 ) 


The last integral is a source term integrated over the enclosed area (drdx). 
Evaluating the flux balance Equation (10) on an arbitrary flow element yields 


- (0%);- p\) E v(\)„v 

k = 0,1, 2 , . . . NH 

where the r term has been approximated by its value of the center of earh flux 
cell . 


\l A s - % ** - 0 < 1] ' 

7 s 



Figure 8 - Cartesian Flux Volume Element Description 


Equation (11) is a discrete approximation whose subscripts refer to the 
respective faces of the elemental volume shown on Figure 8a. Central 
differencing of the flux terms and collecting the contributions at each node 
results in a tridiagonal equation system in terms of ci^. 
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(12a) 


where 


A ’ 9 i,j-l,k + B j 9 i , j , k + C j 9 i , j + l,k = W j,k 

k = 0,1,2,. ..NH 


A j 


B j 


S 

&r 


i ( a e * A w) * i. ( a m * A s) 


A 

c. = 3 

* a r 


a, k = - A w ( 3 2 g, j , k - ^ ( 3 ,! g in , k 
j,k — 1 1.J.K AX 1 



(12b) 

K 2 aA 

(12c) 

r ij 

( 12d) 


( 12e) 


Closure to the problem formulation requires application of boundary conditions 
on the boundaries of the domain of integration. Referring to the computational 
outline in Figure 7, the no flow condition will be imposed on all solid 
surfaces. The upstream flows normally are defined in terms of the engine 
discharge conditions, but flow continuity (Kutta condition) at the lobe 
trailing edge necessitates an alternate approach to avoid overspecifying the 
inlet conditions. These boundary conditions will be explained in the following 
sections. 

Surface No-Flow Boundary Condition 

In a flux formulation, the no-flow boundary condition is implemented as a zero 
flux condition on all solid surfaces. The flux balance for an element 
intersecting the three-dimensional lobe surface is not separable along the 
lines previously derived. In order to render the problem separable, the 
boundary conditions must be linearized in the e direction resulting in a 
modified boundary condition with a source-like term applied on a mean 
axisymmetric surface, R m (x). The surface intersecting computational element 
effectively looks like Figure 8b. The flux balance on such an element must be 
modified by an additional source term representing the "surface" flux (g r ) 
in the r direction. 


rdx - g k 3 dx = 0 

aA 

1,2,.. .NH 
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The areas in the integral^ describe only the external portion of the cell on 
"he side and azini'thal (%A) faces. 


By applying the nean radius approximation to the contour, one can also express 
the lobe surface in terns of separable variables, i. e. , 


NH NK 

f(x,e) = R m (x) + E x k (x) h k (e)= E x k (x) h k (o) (14) 


where \ 0 (x) is the ax i symmetric modal shape and may be used as a mean radius. 

A Fourier moment analysis couples the e dependence of the contour with dependent 
variables. The unknown coefficients \ ^ are determined by using the angular 
definition shown on Figure 9 and arbitrarily defining the potential within a 
constant B^. 


X 


k 


(x) 


B 


k = 


1 


Figure 



f(x,e) cos 



do 



Substituting the modal description of the geometry into the no-flow condition 
and taking advantage again of the orthogonality of alternate Fourier modes for 
the g^'s on the boundary, reduces the surface boundary condition to a flux 
along a mean surface, which is given by 


9 r (x, R ) = \ k (x) + l 
k m * ?r 


2 Z 2 Xn ( *> 3j' x -V p M -j (15) 

m j n L J 

k = 0,1.2,. ..NH 
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where 6-jj is the Kronecker delta function. The first term on the right side 
is the primary contribution to a given mode. The second term represents a 
coupling of the different modal solutions due to the product f e tf 0 in the 
no-flow equation. If parameter e£ is small, this term is absent. Observe 
that Equation (15) is a mixed type boundary condition. Therefore, the solution 
algorithm cannot explicitly determine the gk(x,r)'s from the given boundary 
conditions. This problem of coupled modal equations is alleviated by lagging 
the alternate modes in the iterative solution that will be discussed shortly. 

The modal solutions for the functions g|<(x,ri can be discriticized along to 
the cylindrical mesh to yield 


-tfV A w + 0 2g x A E + 9 r A N 
F x kw w H x kE 11 r kN N 


ks " k 
k = 0,1,2, ...Mil 


a s - g,, V x - 'Lh ** 

R 


( 16 ) 


m 


where A denotes the exterior portions of the areas of the intersected cells 
and aE the exterior area of its azimuthal face, see Figure 8b. In general, a 
surface element includes a region of flow and a region interior to the body. 
Along the mean surface however, an element includes a core flow and fan flow 
region. Correct flux balancing in each region is expedited by tracking the 
potentials of each flow separately across the surface. 

Di scri ticizi ng the boundary Equations (15) and (16) is complicated by the 
mixed mode term in the "surface" flux. In order to solve the equations 
numerically, the coupled term is introduced into the Wj f k right side tern 
and lagged in the iterative solver, i. e.. 


“j.k - 

M j,k * V x 


(17a) 

■j = B J 

* \ 12 , G i 

2R^ 


(17b) 


where Gj and G 2 are the coefficients from Equation (15) of the kth mode 
and the mixed modes (minus the kth term) respectively. 


Solution Algorithm 

The governing equations and boundary conditions reduce the analysis problem to 
a system of linear algebraic equations in terms of the gk(x,r)'s and 
correspondingly the velocity potential. These equations are solved iteratively 
using a successive line over-relaxation procedure (SLOR). To optimize the 
calculations, a grid halving algorithm is utilized. In such an algorithm, the 
previous coarse grid solution is interpolated onto the next or finer grid as 
its initial guess. The solution convergence on each grid is monitored by 
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tracking either the residual, defined as the nornalized error in the nass 
conservation equation at the nth iterate, or the jump in potential at the 
splitter trailing edge. This latter variable is the "flap" vorticity discussed 
in the introductory renarks. Typically the residual is a monotonical ly 
decreasing function that, in a SLOR schene, achieves only a two orders of 
magnitude reduction per grid. Finally, the results are displayed in terns of 
the pressure coefficient defined for each stream relative to its own upstream 
dynamic head. 

Inlet Flow Boundary Conditions (Compound Choked Flow Analysis) 

Flow requirements for cruise engine operating conditions are typically 
determined by engine power settings as well as by inlet and nozzle exit areas. 
Analysis of the localized mixer/nozzle problem requires a completely specified 
set of boundary conditions; however, precautions must be taken to avoid 
specifying inconsistent fan and engine core flow conditions such that the 
Kutta condition at the lobe trailing edge will not be maintained. For example, 
if the onset flows each were specified in terms of p 0 , T 0 and ;1, the 
problem could be over-defined. Use could be made of the nozzle exit flow 
conditions in conjunction with the wake contact surface matching conditions to 
determine uniquely the fan and core flow requirements. The problem is further 
complicated by trying to insure that the analysis would include both choked 
and unchoked conditions at the nozzle exit plane. Table I details duct 
operating conditions for several nixed flow installations. The table indicates 
that the nozzle exit operating condition is dependent on the mission profile 
of the aircraft and that choked and unchoked conditions are indeed possible. 
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TABLE I 

TYPICAL MIXER DUCT OPERATING CONDITIONS 
( Standard Day) 



E 3 

JT8D-: 

■09 (Super SO) 

JT15D 

(no scallon) 


SLTO 

SLTO 

Avg. Cruise 

SLTO 

T amb 

543.7 

519.0 

393.8 

519.0 

Pamb 

14.7 

14.7 

3.46 

14.7 

/ 

Tof 

625.6 

627.5 

560.2 

602.5 

Fan j 

P 0F 

22.45 

24.36 

9.33 

22.04 

' 

Wf 

1157.2 

299.1 

120.1 

55.5 

1 

't 0 p 

1490.0 

1276.0 

1221.0 

1534.0 

\ 

Primary < 

P 0P 

21.65 

25.60 

10.1 

21.09 


[h p 

159.6 

161.4 

70.5 

21.1 

AF 

36.01 

1080.0 

1080.0 

173.06 

Ap 

908.0 

648.3 

648.3 

130.9 

a N 

-- 

1096.7 

1096.7 

304.0 


/ t on 

737.3 

864.1 

814.3 

919.4 


J p 0N 

22.20 

24.76 

9.60 

21.46 

Nozzle 

j 

1 aMI X 

85.0a 

92.0a 

92.0a 

84.0a 


p OM/ p oo 

1.50 

1.68 

2.78 

1.46 

Condi tion 

Unchoked 

Unchok 

?d Choked 

Unchoked 


15 


The analysis developed is based on a compound nozzle model, wherein the 
curvature induced traverse pressure gradients are neglected and the static 
pressure is only a function of axial position. Approximating the nozzle exit 
flow by a quasi-one dimensional flow, the continuity equation for a single flow 

d (pvA(x)) =0 (ig) 

and the isentropic forms of the energy equation can be combined to give 


dA = A 
Cfx 7 



1 \ d In p 

) ^x 


This equation has been generalized 11 for multi stream flows as follows 


( 19 ) 


= N /i 

^ Vtfj 


\ d In p 
) dx 


( 20 ) 


where quasi-one dimensional matching conditions across each dividing 
streamline prevent mixing but allow pressure communication ( p = p) across the 
flow. But area continuity 


therefore, 


*T = 2 A i 
i 


dA. T = £ <*A. 
dx i dx 


v d In p 
dx 


where v , the compound flow indicator is 


Reference 10 also 
to the sign of v , 



1 


) 


demonstrated that the nozzle flow exit state corresponds 
i.e., 


( 21 ) 


v > 0 compound subsonic 

v = 0 compound sonic 

v < 0 compound supersonic 

Equation ( 21 ), in conjunction with the definition of the local mass flows, in 
terms of their stagnation properties. 


w. 


= w. 


<v 


V 


V 


( 22 ) 
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and the isentropic relation are sufficient to provide closure of the problem. 
For a dual stream flow, equation (21) becomes 

„.*1 L -1 * A " L -1 (23) 

n "in ^2 »2n 

The superscripts refer to axial locations defined on Figure 10. This can be 
further reduced using equation (40). In particular, assume that the flow at 
the nozzle plane is subsonic, then 



and the compound flow parameter becomes 



The value of v is not determinate in its present form, since Ai , A 2 are 
not known. An iteration to determine these areas follows. 


(a) Guess value for p^ for given Pq^, A^ n , Aj 1 , AjjjJ. 

(b) Calculate w ^ from following equation given in Reference 11. 



m 

(c) Determine the new value of from equation (22) 

m m 


w 2 ~ w 2 ^ A 2, 17 1, P 02, T 02^ 

but under relax update according to the following 

m _m m _m 

?l = Pi + 0 ( Pi " Pi ) 


(25) 


and continue iteration until 


a P™ within given tolerance. 


(d) Determine 


«1 = Wj 


(a T.p™.Pol. t oi> 
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(e) Determine Aj using equation (22), i.e., solve 


w i = W 1 ( A 1 P„ Pn T n ) for A 


' 1 , 


V °1 


'1 


(f) Calculate v according to equation (21) . 

If v<o , then the solution is complete and the flows are defined in terms of 
the given upstreams data and the flows determined in iteration (a-f). If v=o, 
then the flow is compound choked and p g is not p^. 

For compound choked flow, a dual iteration is performed to find the exit 
pressure for a convergent nozzle, and the splitter exit plane static pressure. 
In this case, it is known that the compound flow indicator v = 0. An initial 
guess for the primary and secondary nozzle flow areas is given by: 

Af = A ; ( Al s /A s t ) 

A 2 = A J ( a|/a|) 

Using equation (24) with p^ replaced by p e , an iteration is performed 
until a value of p e is obtained which gives v = 0. Now the primary and 
secondary flows can be computed from the continuity equation (22) at the exit. 
The secondary flow continuity equation is then used to compute the splitter 
exit pressure, which is then used to compute a primary splitter exit flow. If 
this flow does not match the primary nozzle flow, a new primary and secondary 
nozzle area are computed using the primary splitter exit flow, and the entire 
choked flow calculation is repeated until the flows balance and v = 0. 

Finally, with the flows, pressures and temperatures in both streams known, the 
reference lach numbers at the inlet plane for each stream are computed from 
the continuity equation (22) and the isentropic equations. The velocity gx of 
each strean is therefore directly known. In the present analysis, it is 
assumed that gx = g ;<0 while all higher nodal derivatives at the inlet nlane 
are identically zero. 

Kutta condition - Powered Wake Analysis 


Treatment of the flow downstream of an arbitrary body is complicated by the 
unique interaction of the strean from above and below the body. Even in 
noapowered situations, the lift of flow turning is reflected as a jump in 
potential convected from the sharp trailing edge. This jump will remain 
constant and will follow the trailing stagnation streamline. In the situation 
of power addition, streams are assumed locally irrotational with different 
energy levels. They inviscidly interact through the local potential jump that 
is determined from the basic consistency conditions across a "contact 
discontinuity," i.e., static pressure match on the wake or vortex sheet Sw and 
streamline slope continuity 

7 . T = o on S v/ , S w 
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Classical linear theory assumes the wake lies along a constant radius surface 
(Rm) from the trailing edge, and the streamline slope condition is relaxed. 
Flow is permitted to cross through the constant radius. Consistent with the 
surface boundary condition formulation, the ’wake can be modelled by a mean or 
axi symmetric surface that varies with axial position. 



mean radius approx, 
wake streamline 


Consider the mass flux balance on an element that includes an arbitrarily 
oriented mean wake streamline (Figure 8c). Applying the flux balance (16) to 
the upper and lower portions of the (ij) element results in equations that 
assume a mass flux can exist across the given wake contour, thereby allowing a 
lagging of the wake path in the iteration cycle. Adding the conservation 
components together will produce an_equation for the ij tn element where the 
flux contributions across the wake g r u »*- identically cancel each other 
out. If the potentials on both sides of the wake are defined in terms of a 
mean potential and a wake jump [g^] as follows 


9 ijk = \ (g ijk + 9 ijk ) 

[g k ] fj = 


(26a) 

(26b) 


the flux balance for a wake element becomes 


A j 9 i J-l,k + (B j + B j } 9 i,j,k + C j 9 ij + i,k = W j,k + W j,k ' Cg k ] ij (B J “ B j> 


T 


k = 0,1,2, 


NH 


(26c) 


This Equation corresponds to one row in the general (JMxJM) tri diagonal matrix 
for the ith line. 


A* B.\ C.\ 

j-1 j-1 j-1 


a j c j u 

a j 


b j 


c j 


(27a) 
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The structure of the matrix is equivalent to the following (JM+1) x (JM+1) 
tri diagonal matrix 


! 






0 







(27b) 


The two new (ijth) equations are similar to individual wake element flux 
balances with the added proviso that the A terms are zero. This is equivalent 
to assuming that the wake is a solid boundary since all boundary areas are 
assumed by the program logic to be zero. Under such a format, the program can 
automatically treat these constructed equations in the algorithm and the real 
wake equation is arrived at by contracting or adding the appropriate equations. 


The matrix solution for the i^h row needs an algorithm for defining the 
t g*c J and for updating the wake path. The potential jump, [g^] is obtained 
from the constraints of static pressure and streamline slope continuity. 
Expressing the flows relative to the same reference freestream static 
pressure, the pressure match condition reduces to 


C 


pu 





y n 2 

u u 



2 A Pa 


u ri 


"ii 


(23) 


The pressure coefficient Cp can be taken from the general isentropic 
definition or from a formulation consistent with our linearized algorithm and 
thus separate out the axi symmetric component on each mode C 

Pk 


(x,r) = -2 g ( x,r) 

p k x k 

The streamline slope matching conditions, along the general wake contour 


( 2°) 


1 + 0* 1 * 


77T 


(30) 


can also be simplified by assuming that the axial velocity flux contribution 
is small relative to unity, then 




u 

r 




or for each mode 


If the slope 


dr 



the wake follows the constant radius approximation. 


In order to obtain closure (Kutta condition) Hr or the problem an additional 
assumption is needed. Although the global Kutta condition. 


C 


pu 



(31a) 


needs to be maintained individual modes need their own boundary condition. It 
is proposed therefore that each mode satisfy 


g 


u 

x 


k 



(31b) 


and therefore equation (31) is implicitly satisfied. Equations (28b), (26a, b) 
are then combined to determine the potential .jump 


Cg k ] x = E 1 6ok + E 2\ (32i 

Ideally, the jump in potential along the wake is obtained by integrating out 
axially from the trailing edge along the mean radius as follows 

[9 k ] ij ■ E 1 «ok ,x - x te' * E 2 9 ki * l [9 k ] - E 2 9 k} TE 133 

In the entire preceding discussion, the slope was assumed given. In actuality, 
the path is evolved as the program interates. A simple streamline tracing 
procedure could be used to periodically update the wake path. Another 
attractive approach would be to drive the slope using the nonconservation over 
the local flux cel 1 . 

Equation (1) can now be re-expressed in terms of the separated variables as 
follows 


NH/2 

r= 4S (g 2k AUoo - [g 2k ] Uc*) TE (34) 


where only the odd nodes 


contri bute 


to the net 


induced circulation field. 
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C. Geometry Definition 


General lobe contour definition can be a complex problem even if scalloping 
and scarf angle cutouts are not included. Current commercial designs fall 
within three general categories: radial sidewalls, parallel sidewalls, 

circular arc sidewalls with the remaining segments of the lobe defined*in 
terms of tangentially intersecting circular arc segments. Considering the 
periodic nature of the analysis formulated in the previous section, one must 
limit the geometry capability to axi syimetric duct wall configurations while 
the lobe cross-section must be limited to radial sidewall geometries to avoid 
multiple valued structures. Using this approach, a mixer centerbody and fan 
cowl can be defined and replicated using the BYU Movie Three-Dimensional 
hidden line graphics program, to produce Figure 10. 



Figure 11 Three Dimensional Display of JT8D-209 Mixer Geometry 


A measure of the feasibility of the Fourier decomposition method can be made 
by considering a finite number of harmonics and comparing the reconstructed 
lobe versus the given lobe definition. A equal angle series construction 
permits use of fast Fourier scheme (FFT) to evaluate the Fourier coefficients 
The reconstructed lobes are shown in Figure 12 for MH = 1 , 3 and 9. The base 
contour was generated analytically using 50 points. The modal analysis, 
however, does not demonstrate the perfect agreement. The largest excursions 
occur at points where the curvature changes instantaneously. Whereas this 
aspect is commonly found in square wave reconstructions, where a lame number 
of nodes must be used to obtain an accurate wave representation, the si ight 
lobe contour mismatch should not be significant in establishing the "flap" 
vorticity field. The sensitivity of the analysis to the number of modes, 
however, will be shown below to be a critical factor in determining the’nean 
radius, the leading tern obtained from the FFT analysis. 
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SECTION 3.0 

RESULTS AND DISCUSSION 


IivJr? 1 TnK S ^ F - 0MIX) d ®! cn ’ bed ™ the previous section has been applied to 
several lobed mixer configuration. When Xu = 0 (k=l, 2....NH) the zeroth 

mode solution will describe the powered flow over a completely axisymmetric 
configuration. An initial calculation of a olanar mixer lobe in a straioht 

bv C ^nn^H^- ented t0 . call ' bra te the analysis. Planar conditions are simulate* 
by considering an axisymmetric geometry at large radius. In such a situation 
comparison calculations can be made with an available analytically constructed 

dniiMp?* t i at S1rnulates an isolated planar mixer lobe using distributed 
doublets along a mean planar surface** (PLANMIX). Figure 13 shows a sideview 
of the configuration analyzed. The duct walls are defined sufficiently far 
fr« the l°be S (H = 8") so that any interactions would he Si JimeHcal 

thl^P w^i? nS d ® mo J s , brate * hat there 1S no potential interaction effect due to 
these walls. The lobe surface in both calculation methods is generated from a 

fwL?H $ T h l he A ref ° re 3 single moda1 (NH = 1} sol *Jtion models the 
tlowfield. In this study no power addition effects 

are presented for the reconstructed (as a function 
perturbation velocity on the lobe surface. Figures 
comparison calculations between FLOMIX and PLANMIX 
several azimuthal cuts running from the lobe crest 
the lobe trough ( 0* = 1.0). 
its trailing edge. The nrol 

approxi nations at the trailing edge, i.e., the Kutta condition is satisfied 
h lnl r 8 velocity perturbation goes to zero. Althounh both methods arn 
ba^ed on linear theories, slight differences should be expected since the 

an 1 nverse method evaluating singular integrals nune-icall v 
as input v/mle the present method is a finite difference flux volume scheme. 


are considered. Predictions 
of e) components of the 
14(a), (b), (c) show 
the three components for 
- — ( 0,= 0) to the inside of 
T.ie axial scale runs from the lobe leading edge to 
'*iles reflect the effect of the linear theory " " 



igure 13 Sideview Representation of Planar Mixer/Wind Tunnel Geometry 
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1 .05 



Figure 14a Comparison Calculations of Axial Velocity for a 
Symmetric Planar Mixer in a Planar Duct 



Figure 14b Comparison Calculations of Vertical Velocity for a 
Symmetric Planar Mixer in a Planar Duct 



Figure 14c Comparison Calculations of Spanwise Velocity for a 
Symmetric Planar Mixer ir a Planar Duct 
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More realistic applications of the method can be found by considering flight 
type configurations. In particular, two specific powered applications are 
presented due to the "benchmark" nature of their experiments: (1) the Energy 
Efficient Engine (E 3 ) configuration 29, which is an 18 lobe forced mixer^,Io j 
and (2) a JT8D-209 12 lobe forced mixer^. The E^ configuration, shown in 
Figure 13, is well suited to the present formulation in that although it was 
modeled for a modern high bypass engine, it was designed specifically for code 
verification, i.e., no scalloping or scarfing of the lobes was used and 
extensive surface static pressure surveys were made in the lobe region of the 
mixer. In contrast, the JT8D-209 forced mixer is a higher penetration 12° 
scarf angle design typical of first generation low bypass applications. While 
this geometry is not strictly suitable for code comparison, the experimental 
data included LDV profile measurements of all components of velocity at the 
lobe trailing edge. Satisfactory modeling of these lobe cross-sections, even 
for the JT8D-209 high penetration lobe (Figure 12) is possible with only ten 
Fourier terms, however, a more definitive approach to defining NH is given 
below. 



TE VIEW SIDE VIEW 


Figure IF Schentic Representation of E 3 Configuration ?9 Lobed lixer 
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Postprocessing of the solution produces surface Mach number definition and 
quantifies the level of "flap" vorticity or circulation through the potential 
jump at the lobe trailing edge. The surface solutions, although strictly 
determined along the approximate mean surfaces can be interpreted, to first 
order, as solutions on the actual surfaces. By examining the potential jump 
variable, an effort will be made to separately quantify the contrihutive 
effects of mean flow turning, power addition and lobe aspect ratio 
(penetration angle) on the overall mixing process. 

Numerical solutions initially have been made for the E 3 configuration 29 
lobed mixer and comparisons have been made with test data measured on a full 
scale model at FluiDyne Engineering Corp. The experimental cruise flow 
conditions were characterized by power settings* of aP 0 /P 0 s = 9.094 and 
aTq/Tcs = l- 3n . Computational simulation of the measured flow conditions 
is obtained by setting the upstream total conditions in each stream as well as 
settings the exit static pressure. The quasi one-dimensional choked flow 
analysis correctly sets up a choked flow at the nozzle exit plane and sets the 
inlet flows and Mach numbers to within 10% of those measured by the facility 
flowmeters. Solution accuracy, relative to the level of modal approximation, 
however, is to be answered. 

Table I summarizes the individual modal potential jumps [g|<] obtained from a 
series of calculations in which different levels of modal approximation (NH = 
0,1, 2, 3...) were used to simulate the E 3 mixer. The calculations were 
initially made for an axisymmetric configuration where the effect of flow 
turning of the mean radius (-0.093) and power addition (-0.012) could be 
identified. With the mean radius having a positive trailing edge angle, these 
results demonstrate that a positive circulation corresponds to a clockwise 
rotation. The tabulated calculations varied both the number of harmonic terms 
(NH) used to represent the © dependence of the velocity potential and the 
number of Fourier terms (NF) used to represent the e dependence of the lobe 
surface. The tabulated results indicate that while each approximation produces 
different total and modal potential jumps, the individual modes approach fixed 
values as more terms are included. Although the higher modal solutions 
converge much more rapidly than the leading terms, the additional storage 
required for these higher order terms is impractical. A study of the governing 
differential equations indicates that the modal couplino is extremely weak Tor 
the axisymmetric mode and that its primary driving term is the mean radius, 
determined from the Fourier analysis. Since higher order solutions (MH large) 
contribute little to the total circulation field, one can neglect these 
equations while retaining the addition terms for an improved Fourier (NF) 
definition. For example, calculations with NH = 3, NF = 5 yield the same modal 
jumps as MH = 5, MF = 5. In the comparisons described below the solution 
parameters NH = 5, NF = 18 were used as representative of an "asymptotic 
solution. 


* The subscript s refers to reference conditions in the secondary or fan 
stream. 
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TABLE I. POTENTIAL JUMP COMPARISONS [g k ] 


Power Total 

NH MF Addition Jump k=0 


0 axi (18) 

no 

-0.093 

-0.093 

0 axi (18) 

Yes 

-0.118 

-0.118 

1 1 

Yes 

-1.396 

-1.328 

3 3 

Yes 

-0.527 

-0.467 

5 5 

Yes 

-0.353 

-0.292 

3 5 

Yes 

-0.367 

-0.293 

5 9 

Yes 

-0.263 

-0.187 

Asymptotic Value 

-0.191 

-0.118 


(MH=5, NF=18) 


1 

2 

3 

4 5 

-0.068 

-0.084 

0.018 

0.006 


-0.C87 

0.013 

0.007 

0.006 

-0.087 

-0.088 

0.013 

0.008 

0.009 

-0.004 -0.003 

-0.089 

0.008 

0.010 

-0.002 -".003 


MACH 

NUMBER 



AXIAL DISTANCE ~x/L 


Figure 16 Surface Mach Number Comparisons for E 3 Duct Walls at 
Lobe Crest Orientation 
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Figure 17 Surface Mach Number Comparisons for E 3 Mixer Lobe at 
Lobe Crest Orientation 


MACH 

NUMBER 



Figure 18 Surface Mach Number Comparisons for E 3 Duct Walls at 
Lobe Trough Orientations 
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MACH 

NUMBER 



Figure 19 Surface Mach Number Comparisons for Mixer Lobe at 
Lobe Trough Orientation 


Surface lach number calculations, compared to measured data within the lobe 
region, are shown on Figs. 15-19. Figures 15, 17 present comparisons made for 
an azimutoal cut aligned to the lobe crest while Figs. 18, 19 show comparisons 
for the lobe trough orientation. Initial calculations indicated that the 
quasi -one- dimensional boundary condition set the flow and inlet conditions 
approximately 10 percent too high, therefore, a nozzle flow coefficient (C„ 

- 0.94) was introduced to adjust the inlet flow conditions. With this / 
modification, the figures show substantial agreement between analysis and 
data. The axisymmetric fan nozzle and centerbody comparisons shown on Figures 
and 18 show little angular variation and are largely one dimensional in 
behavior. The lobe surface solutions, however, show substantial e dependence. 

A major discrepancy is noted near X/L = 0.50, the cross-over point for the fan 
and core flows. A Mach number pulse at this point is produced primarily by the 
first modal solution in a manner similar to Figure 12a, observed in the 
previously discussed planar mixer case. Viscous interaction effects should 
decrease the analytically predicted gradients and reduce this mismatch. 

The FLOMIX analysis has also been applied to the JT8D-209 lobed mixer 
configuration. This configuration was previously studied at IJTRC in a scaled 
model test but still simulating the hot flow (APq/P 0S = 0.044, AT 0 /T 0S 

= 1.617) full scale engine cruise conditions. Analysis of this conf iguration 
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is complicated by the additional effect of a 12° scarf angle. While no surface 
or flowfield details were measured within the lobe region of the duct, LDV 
measured velocity components were obtained at the lobe trailing edge and in 
the mixing duct. Any analysis comparison with data, however, requires an 
interpretation of the small disturbance solution off the mean radius. 
Examination of the Mach number indicates that axial component is largely 
one- dimensional and has little variation with radius, within each stream. 
Displacing the flow relative to the physical lobe trailing edge results in 
comparisons for the radially measured axial component of Mach number, shown on 
Figure 20 at the crest location and on Figure 21 for the trough location, both 
figures show substantial agreement between analysis (NH = 5, MF = 9) and 
experiment. , 



NORMALIZED RAD AL DISTANCE 

Fiqure 20 Comparison of Axial Mach Number at Lobe Trailing Edge 
Plane for JT9D-209 Mixer, Crest Orientation 


d'= i.o 



NORMALIZED RADIAL DISTANCE 


Figure 21 Comparison of Axial Mach Number at Lobe Trailing Edge 
Plane for JT8D-209 Mixer, Trough Orientation 
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A fuller interpretation of the flowfield at the lobe exit plane must be 
inferred from slender body theory 14,15^ whereby the outer potential is 
determined as a function of x from the solution of an axisymmetric problem and 
the inner potential is determined as a solution of the two-dimensional 
Laplace s equation in the cross plane (r,e). This philosophy will be used with 
the modal axisymmetric solutions to infer the cross flow or inner solution by 
viewing the potential jump as the doublet source equivalent obtained from the 
outer solution. 
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A. Geometry Definition 


SECTION 4.0 

ANALYSIS IMPLEMENTATION 


A generalized procedure to obtain the coordinates of an arbitrary forced mixer 
application can be extremely complex. Although an engine centerbody is usually 
an axisymmetric surface, the outer cowl, which is initially also axisymmetric , 
can transition to a high AR rectangular cross section at the nozzle exit plane 
(Figure 22). Commercial applications are typically fully axisymmetric, but 
military applications can include such ducts. Superellipsoidal coordinates can 
be used to analytically approximate such circular to rectangular transition 
ducts. 



Figure 22 Engine Exhaust Transition Duct. 

Lobe contour definition is however a more complex nroblen. Even if scalloping 
and scarf angle cutouts are not modelled, lobe contours can still be very 
general. Current commercial designs fall within three general categories; 
radial sidewalls, parallel sidewalls, circular are sidewalls (See Figure 23). 
The remaining segments of the lobe are then defined in terns of tangential Iv 
intersecting circular a**c segnents. 
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Considering the periodic nature of the analysis fornulated in the previous 
section, one must limit the geometry capability to axisymmetric duct wall 
configurations. The lobe cross-section must also be limited to radial sidewall 
geometries to avoid multiple valued structures. Lobe coordinates are 
analytically evaluated by dividing the lobe into three segments (shown on 
Figure 24) . 



Figure 24 Analytical Breakup of Radial Sidewall Lobe 

Using this approach, the baseline JT8D-209 model 12 mixer is defined interms 
of the radii given on Table II, where the appropriate radii and the sectional 
locations are shown on Figure 25. Finally, after generating the mixer 
coordinates, one can replicate the lobe and display, using the BYU Movie 
Three-Dimensional hidden line graphics program, the lobe portion (Figure 26) 
as well as the complete forced mixer configuration (Figures 27, 28). 


Table II - JT8D-209/12C Lobe Mixer Definition 


AXIAL SECTION 

A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

M 

N 

X 

.655 

.951 

1.236 

1.593 

1.837 

2.093 

2.284 

2.713 

3.141 

3.570 

3.998 

4.306 

4 .613 

5.399 

R r i n V » t Icy 

2.729 

2.719 

2.694 

2.642 

2.593 

2.517 

2.449 

2.276 

2.103 

1.929 

1.756 

1.632 

1.508 

1.190 

Rcr OW n 

2.503 

2.503 

2.503 

2.503 

2.526 

2.593 

2 . 662 

2.831 

2.994 

3.105 

3.160 

3.167 

3.146 

2.997 

D 








.301 

.313 

.310 

.296 

.279 

.258 


E 

Rrv Tangency Radius 


2.512 

2.295 

2.088 

1.891 

1.757 

1.628 


F 

Rhc.n 


.313 

.313 

.313 

.313 

.313 

.313 


G 

Rile • n 

- [Rrv + Rc 

r o w n J / 2 




2.243 

2.243 

2.243 

2.243 

2.243 

2.243 


H 

Re r own 



.299 

.274 

.250 

.235 

.233 

.239 

J 

R Crown Tangency Radius 


2.56 

2.760 

2.902 

2.975 

2.904 

2.957 


K 

Rrv Origin 

Radius 




.963 

. 354 

.303 

.247 

.207 

. 197 

.173 


L 

R Crown Origin Radius 

. 300 

. 300 

.278 

.255 

.241 

.239 

.244 
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Figure 26 3-Dimensional Movie Plot of Input-Generated JT8D-209 Mixer 




Figure 27 3-Dinensional Movie Plot of 
Complete JTQD-209 Mixer 
Geometry 


Figure 23 3-Dimensional Movie Plot of 
Complete JT3D-209 Mixer 
Geometry (Looking Upstream) 


Finite Fourier Series Decomposition 


Lobe contours are constructed using tangential intersection of radial lines 
with a variety of circular arc segments. The contours are assumed symmetric 
about the either the crest or trough of the lobe. The half lobe is then 
subdivided into NINT equal angle segments. At these points (Rf.e-j) the 
lobe can be represented by a finite Fourier series, where the number of 
Fourier Terms is determined by the (NINT+1) constraints. A measure of the 
feasibility of the Fourier decomposition method formulated in Section II can 
be made by considering a finite number of harmonics (NH = NINT/2) and 
comparing the reconstructed lobe versus the lobe definition array from the 
given input data. 


The definition of an arbitrary point on the lobe surface is given by the 
following series 


where 


F = R.(e.) = j 1 A. cos (ka.) + B|< sin (kai) 
J J k=o [_ J J 


+ R , 


CM 


2*0. / 2e o = ir0j/o o 


0 < ©j < 2 o q 



Evaluation of the unknown coefficients is possible by forming 2NH moment 
relations that take advantage of the orthogonality of the basic Fourier 
functions, 
e.g. 


2tt 


S t 


AF cos( ia) da = 


NH 

2 

k=l 


2* 2ir 

^A^ cos ( ka) cos( ia) da + ^/^k Sln *a) cos( id) da 


i = 0,1,2, NH 


3 ' 


Using trigonometric integral identities to 
obtains 


B 

A 




F sin ( k a) da 
F cos (ka) da 


simplify the above equation one 
k = 0,1,2, ...NH 
k = 0,1,2, ...NH 


An equal angle series construction permits use of fast Fourier scheme (FFT) to 
evaluate the unknown coefficients. A more reasonable contour definition, 
however, can be obtained by using an equal arc length definition of the input 
data. In such a formulation, the coefficients must be evaluated using a more 
conventional matrix inversion procedure. 

The truncated Fourier series has been applied to the three lobe families shown 
on Figure 23. The limit on angular double valuedness for the parallel and 
curved side wall geometries was eliminated, only for this geometry study, by 
using an equal arc length decomposition of the lobe surface. The reconstructed 
lobes are shown on Figures 29, 30, 31 for NH = 3 and 5. The base or reference 
contour was generated analytically using 50 points. One can see that 
relatively few modes are required to produce a good surface representation of 
the parallel and curved side wall cases. The radial analysis however does not 
demonstrate the same degree of agreement. The largest excursion occurs at 
points where the curvature changes instantaneously. This aspect is commonly 
found in square wave reconstructions, where a large number of modes must be 
used to obtain an accurate wave representation. The relative mismatch is, 
however, enhanced by the magnification of the abcissa scale. 

The level of error should be compared however not only to the analytic contour 
reference but also to the manufacturing tolerance of the mixer. Mixer lobes 
can be accurately fabricated by three-dimensional machining of a solid mass, 
however, this is prohibitively expensive. Our JT8D-200 mixer experience, also 
used for the JT8D benchmark mixer (Ref. 5), shows that a typical mixer can 
have dimensional variations as large as 0.4 inch. For example, a calibration 
of two adjacent lobes on the benchmark mixer (Figures 32, 33) shows 
substantial lobe to variation. It therefore seems reasonable that only 3-5 
modes will be necessary to provide adequate three-dimensional contour and flow 
reconstruction. In production practice the lobe-to-lobe variations seem to 
compensate as long as lobe flow area is maintained. 
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Figure 30 Comparison of Fourier Modal Reconstruction of Lobe 
Contour for Radial Sidewall Model of JT8D Lobe 
at STA 1 3,X/L= 93% ( X Fourier, o Analytic) 
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Figure 31 - Comparison of Fourier Modal Reconstruction of Lobe 

Contour for Radial Sidewall Model of JT8D Lobe 
at STA 7,X/L= 45% (...X Fourier, o Analytic) 






Figure 33 Circumferential Profile of JT8D-209 baseline Mixer 
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B. Linking Diagram for FLOMIX Program 


MAIN — ! 


ECHO 


— FCRMIX 


j- RWINPT 
— PLOT 

- IOBGEN 
— EQAMGS 
— FOREGO 

- FORVAL 
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- PLOT 
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PLOT 
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GRID 
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— eCSETS 
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IHITL 
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■- TRIOI 
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MATRIX 


— UPDATE 
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c. Input Description for FLOMIX Program 

The format of each input item is identified by type. 


Type 


A 

F 


I 


Example 
TEST CASE 


-23.64 


-.2364 E2 
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Description 

Alphanumeric: Any keyboard characters are 

specified within the given field. 

Floating Point: Decimal fractions including 

decimal point are specified anywhere within 
the given field. Positive values are assumed 
unless the value is preceded by a minus sign. 

Scientific: notation may be substituted by 
specifying a decimal mantissa (as above) and 
a right adjusted base 10 multiplier preceded 
by a symbol E. 

Integer: Right adjusted whole numbers (no 

decimal point) are specified within the 
given field. 


Blank F or I fields are set equal to zero, but blank A fields are set equal to 
blank characters. 


Card Type 1 

Column Item Type Description 

1_52 AXITLT A Title for plotting routine 


Card Type 2 

Column Item Type Description 

1_8Q TITLE A Title of case for printout 


Card Type 3 


Column 

Item 

Type 

Description 

1-5 

NLOBE 

I 

No. of mixer lobes 

6-10 

fJHARM 

I 

Mo. of harmonic solutions 
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Card Type 4 


Col umn 

Item 

Type 

1-10 

PT1 

F 

11-20 

TT1 

F 

21-30 

GAM1 

F 

31-40 

Rl 

F 

41-50 

P1NF 

F 

Card Type 5 

Column 

Item 

Type 

1-10 

PT2 

F 

11-20 

TT2 

F 

21-30 

GAM2 

F 

31-40 

R2 

F 

Card Type 6 

Col umn 

Item 

Type 

1-5 

IMXFIN 

I 

mesh 

6-10 

JMXFIN 

I 

mesh 

11-20 

XMIN 

F 

21-30 

XMAX 

F 

31-40 

RMIN 

F 

41-50 

RMAX 

F 


Description 

Total pressure of primary flow, psi 

Total temperature of primary flow, °R 

specific heat ratio of primary flow 

Rl, 9 as constant of primary flow 
ft=lb/lbm-‘'R 

Static pressure of external or ambient 
flow, psi 

Description 

Total pressure of fan flow, psi 

Total temperature of fan flow, "R 

7 2 , specific heat ratio of fan flow 

R2, gas constant of fan flow, 
ft-lb/lbm-' , R 

Description 

No. of axial (I) grid lines on finest 

No. of radial (J) grid lines on finest 

X m in» minimum axial grid location 
X max> maximum axial grid location 
*Rmin» minimum radial grid location 
* R max» maximum radial grid location 


♦Coordinate data for centerbody and fan cowl geometries must be specified 
within these radial grid limits. 
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Card Type 7 

Maximum of five grid halvings are possible, with last being on finest mesh. 
Setting NGRIDS=1 will calculate only on finest mesh. 


Col umn 

Item 

Type 

Description 

1-5 

NGRIDS 

I 

No. of grids used in mesh halving 

6-10 

MAXSWP(l) 

I 

Maximum No. sweeps on grid 1 without 
satisfying its convergence tol . 

11-15 

MAXSWP( 2) 

I 

Maximum No. sweeps on grid 2 

16-20 

MAXSWP( 3) 

I 

Maximum No. sweeps on grid 3 

21-25 

MAXSWP(4) 

I 

Maximum No. sweeps on grid 4 

26-30 

MAXSWP 

(NGRIDS) 

I 

Maximum No. sweeps on grid 5 


Card Type 8 


Column 

Item 

Type 

Description 

1-10 

RLXSUB 

F 

Over Relaxation factor for SLOR 
(Recommended=l .7) 

11-20 

PHIXT 

F 

Factor controlling 0xt artificial 
time term 

21-30 

PHIR 

F 

Factor controlling 0^ artificial time 
term 

Card Type 9 




Extrapolated Relaxation 

Parameters 

(Set to zero for no extrapolation) 

Col umn 

Item 

Type 

Description 

FREE FORMAT 

INSTEP 

I 

Minimum No. iteration sweeps before 
forced extrapolation ( defaul t=250) 


INDJMP 

I 

Minimum No. iteration sweeps between 
extrapolation cycles ( defaul t=25) 


TOLJMP 

F 

Tolerance factor (defaul t=l. 5) 


NSWAVE 

I 

No. of sweeps used in calculating 
average residual for extrapolation 
( defaul t=5) 
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Card Type 10 


Output control 

parameters 



Col umn 

Item 

Type 

Description 

1-5 

I DUMP 

I 

=0, Minimum output (Input, 
Fourier analysis, surface Cp) 
=1, + Surface Mach No. 

=2, + surface g's 

=3, + diagnostic output dump 

6-10 

MTH 

I 

Number of azimuthal output planes 
at lobe trailing edge 

11-15 

NR 

I 

Number of radial output points at 
lobe trailing edge 

16-20 

IPLOT 

I 

=0, No output plots 
=1, Summary plots 
=2, + Final grid Cp solution 
=3, + potential and Cp for all 
gri ds 


Card Type 11 

Specify NSTA cards of input data (Program will read to end of file) 


Col umn 

Item 

Type 

1-10 

STA(I) 

I 

11-20 

RID(I) 

I 

21-30 

REL(I) 

I 

31-40 

ROD(I) 

I 

41-50 

RFV ( I ) 

I 

51-60 

XPCNT(I) 

I 


Description 

Axial Station, in 

RID, Centerbody/plug radius, in. 

R crown, lobe crown or maximum radius, 
in. 

ROD, Fan cowl radius, in. 

RFV, lobe fan valley or minimum radius, 

in. 

a/e , Percent of half lobe angle to 
radial side-wall plane (input some 
dummy value for axi symmetric problems) 
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0. Sanple Input (JT8D-209 Forced Mixer) 


INPUT ECHO (CARO IMAGE ) 

10 20 30 *0 50 60 70 00 

COL 123656759 123656769 123656769 123656789 123656769 123656789 123656789 123656789 



1 

FLOMIX 

TEST CASE 






2 

JT6D MIXER 

NOZZLE SMOOTHED DATA (66 EQUALLY SPACED 

POINTS ) 


3 

12 6 







6 

38.1 

1360.0 

1.6 

53.3 

16.5 


CARO 

5 

36.5 

519.7 

1.6 

53.3 

1.0 



6 

129 129 

0.0 

10. 0 


6.3 



7 

3 500 

600 300 






6 

1.50 







9 

o e 

0 0.0 





CARO 

10 

1 9 

20 






11 

0.30000 

1.37000 

2.50306 

3.59961 

2.72635 

.5 


12 

0.50000 

1.37000 

2.50220 

3.60273 

2.73196 

.5 


13 

0.70000 

1.37000 

2.69995 

3.59808 

2.73019 

.5 


16 

0.90000 

1.37000 

2.69682 

3.56596 

2.72256 

.5 

CAPD 

15 

1. 10000 

1.37000 

2.69372 

3.56682 

2.70907 

.5 


16 

1.30000 

1.37000 

2.69210 

3.56987 

2.68888 

.5 


17 

1.50000 

1.37000 

2.69592 

3.53216 

2.66106 

.5 


16 

1.70000 

1.37000 

2.50770 

3.51697 

2.62676 

.5 


19 

1.90000 

1.37000 

2.56077 

3.50659 

2.57560 

.5 

CARO 

20 

2.10000 

1.37000 

2.60666 

3.69627 

2.51038 

.5 


21 

2.30000 

1.37000 

2.66690 

3.68539 

2.63535 

.5 


22 

2.50000 

1.37000 

2.76950 

3.67558 

2.35800 

.5 


23 

2.70000 

1.37000 

2.86637 

3.66676 

2.27979 

.5 


26 

2.90000 

1.37000 

2.91311 

3.65260 

2.20069 

.5 

CARO 

25 

3.10000 

1.37000 

2.97525 

3.63867 

2.12076 

.5 


26 

3.30000 

1.37000 

3.03066 

3.62366 

2.06016 

.5 


27 

3.50000 

1.37000 

3.07763 

3.60716 

1.95911 

.5 


28 

3.70000 

1.37001 

3.11696 

3.39006 

1.87790 

.5 


29 

3.90000 

1.37027 

3.16233 

3.37185 

1.79665 

.5 

CARO 

30 

6.10000 

1.36066 

3.15591 

3.35311 

1.71566 

.5 


31 

6.30000 

1.33175 

3.16393 

3.33630 

1.63660 

.5 


32 

6.50000 

1.29293 

3.15617 

3.31535 

1.55361 

.5 


33 

6.61300 

1.26135 

3.16600 

3.30670 

1.50800 

.5 


36 

6.90000 

1.157*8 

0.0 

3.27836 

0.0 

.5 

CARO 

35 

5.10000 

1.08316 

o.e 

3.26075 

0.0 

.5 


36 

5.30000 

1.00883 

0.0 

3.26317 

0.0 

.5 


37 

5.50000 

0.93637 

0.0 

3.22710 

0.0 

.5 


38 

5.70000 

0.85979 

0.0 

3.21132 

0.0 

.5 


39 

5.90000 

0.78520 

0.0 

3.19571 

0.0 

.5 

CARO 

60 

6.10000 

0.71067 

0.0 

3.16172 

0.0 

.5 


61 

6.30000 

6.63656 

0.0 

3.16776 

0.0 

.5 


62 

6.50000 

0.56219 

0.0 

3. 15639 

0.0 

.5 


63 

6.70000 

6.68860 

0.0 

3.16216 

0.0 

.5 


66 

6.90000 

0.61623 

0.0 

3.12992 

0.0 

.5 

CARO 

65 

7.10000 

6.36666 

0.0 

3.11636 

0.0 

.5 


66 

7.30000 

0.27501 

0.0 

3.10762 

0.0 

.5 


67 

7.50000 

0.20665 

0.0 

3.09663 

0.0 

.5 


66 

7.70000 

0.13968 

0.0 

3.08636 

0.0 

.5 


69 

7.90000 

0.07633 

0.0 

3.07611 

0.0 

.5 

CARO 

50 

8.10000 

0.01003 

0.0 

3.06597 

0.0 

.5 


COL 123656789 123656789 123656789 123656789 123656789 123656789 123656789 123656769 
10 20 30 60 50 60 70 80 



E. Subroutine Definition for FLOIIX Program 


^WwwwWW 


c 

c 

C WELCOME TO FLOMIX, THE AXI3YMMETRIC POTENTIAL FLOW ANALYSIS. THE 
C GUOERLEY-VON K ARMAN TRANSONIC SMALL DISTURBANCE EQUATION IS THE 
C 60VERNING PARTIAL DIFFERENTIAL EQUATION. THE DISCRETIZED 
C FORMULATION IS BASED ON A FINITE VOLUME CONSERVATION OF MASS 
C FLUXES. THE BOUNDARY COWJITIONS ARE APPLIED ON THE EXACT SURFACE. 

C A NONCONSERVATIVE S«OCK POINT OPERATE is AVAILABLE IN LIEU HF 
C THE DEFAULT CONSERVATIVE OPERATOR. THE CONFIGURATIONS WHICH 
C CAN BE ANALYSED INCLUDE STING-MOUNTED OR FINITE LENGTH NACELLES 
C KITH OR WITHOUT A CENTER BODY. BOTH FREESTREAM OR WIND TUNNEL 
C CONDITIONS CAN BE SIMULATED. AN ACTUATOR DISK MODEL IS USED TO 
C SIMULATE POWER ADDITION TO THE FLOW. GRID REFINEMENT AND 
C EXTRAPOLATION TECHNIQUES ARE USED TO ACCELERATE CONVERGENCE. 

C 

C THE MAIN ROUTINE CONTAINS THE BASIC STRUCTURE OF THE PROGRAM. 

C 

C ROUTINES CALLED: 

C ECHO - ECHOS CARO IMAGE INPUT 

C FORMIX - READS THE SPECIFIED DATA FILE, OR HELPS THE USER CREATE 
C A NEW DATA FILE WHEN ONLY THE CONTOUR IS GIVEN. 

C EBC - CALCULATES THE POSITION OF THE DOWNSTREAM BOUNDARY AND THE 
C DOWNSTREAM MACH NUMBERS FOR EACH STREAM. 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GRID - GENERATES A CARTESIAN MESH FOR THE OENSEST GRID. 

XRSUBS - GENERATES COARSE GRID SUBSETS OF THE DENSEST GRID, 

AS WELL AS CALCULATING COARSE GRID VALUES OF VARIABLES 
SUCH AS IMAX, JMAX * IH, AND ITE. 

BCSET3 - GENERATES THE ARRAYS WHICH ARE USED TO ESTABLISH THE 
EXACT SURFACE BOUNDARY CONDITIONS. 

PGUESS - GENERATES AN INITIAL GUESS FOR THE POTENTIAL, PHI. THE 
GUESS IS OBTAINED EITHER FROM A SAVED SOLUTION, THE 
SOLUTION FROM THE PREVIOUS GRID, OR A GIVEN FUNCTIONAL 
FORM. 

INITL - SETS UP INITIALIZATION FOR EXTRAPOLATED RELAXATION 
FARFLD - GENERATES THE FAR FIELD BOUNDARY CONDITIONS. THIS 

SUBROUTINE IS INCLUDED IN THE RELAXATION SWEEPING LOOP TO 
ALLOW FOR PERIODIC UPDATING BASED ON THE CIRCULATION. 

SLOR - PERFORMS A SUCCESSIVE LINE OVER -RELAXATION. FOR 

EACH COLUMN IN THE CARTESIAN MESH, A TRIDIAGONAL MATRIX 
IS SET UP AND SOLVED. THE CORRECTIONS ARE APPLIED WITH 
AN OVER-RELAXATION FACTOR FOR SUBSONIC FLOW, AND AN 
UNDER -RELAXATION FACTOR FOR SUPERSONIC FLOW. WHEN 
THE PROPER CRITERIA ARE MET, AN EXTRAPOLATION IS CALLED. 
UPOATE - UPDATES CORRECTIONS FOR NEXT SWEEP THROUGH RELAXATION LOOP 
RESCU - CALCULATES RESIDUES FOR EACH SWEEP THROUGH RELAXATION LOOP 
OUTPUT - CALLS THE VARIOUS OUTPUT SUBROUTINES THAT APPLY 

TO EACH GRID. THESE INCLUDE THE CP PLOTS, THE MASS 
FLOW INTEGRATIONS, AND THE DRAG CALCULATION. 

FINISH - PROVIDES A SUMMARY OF THE ANALYSIS, AND WRITES 
INFORMATIONS TO A FILE FOR FUTURE CALCULATIONS. 

RECOMB - RECOMBINES EACH HARMONIC SOLUTION INTO ONE SOLUTION 
PJMPLT - PLOTS THE JUMP IN PHI OVER ALL ITERATIONS 

ROUTINES CALLED BY: 


SUBROUTINE ECHO 

c 

C CARD IMAGE INPUT ECHO 
C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
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S^R^IN|_F^MIX< THE, AXITLT ) 

FORMIX DOES THE FOLLOWING: 

U) READS THE LOBE CROSS-SECTION GEOMETRY OATA FROM R. WILEY f S 
PROGRAM. 

(2) FORMS A SERIES OF POINTS AT EQUAL INTERVALS OF ANGLE PHI 
TO DEFINE THE LOBE CROSS-SECTION FOR FOURIER DECOMPOSITION. 

(3) FORMS A FOURIER DECOMPOSITION OF A SELECTED CROSS-SECTION. 

(4) RECREATES THE MIXER CROSS-SECTION USING A SELECTED NUMBER OF 
TERMS OF THE FOURIER EXPANSION. 

(5) GIVES A TEKTRONIX SCOPE PICTURE OF THE MIXER CROSS-SECTION 
AND THE FOURIER REPRESENTATION, 

(6) CALCULATES A TABLE OF THE MIXER CROSS-SECTION VALUES, 
CORRESPONDING FOURIER VALUES, ABSOLUTE ERRORS BETWEEN 
CORRESPONDING VALUES, AND GIVES THE RMS ERROR AND MAX. ABS. 
ERROR FOR THE CROSS-SECTION. 

OUTPUT ARGUMENTS: 

THE - THETA ANGLES 
AXITLT - PLOT TITLE 


ROUTINES CALLED: L06GEN, EQANGS, FOREGO, FORVAL, RMINPT , REPLIC, 

OUT 3D, LAMDAS 

LOBGEN - CALCULATES MIXER LOBE SHAPES AND CROSS SECTIONS 

EQANGS - COMPUTES THE SERIES OF (Z,Y) PAIRS WHICH DESCRIBE THE 

SHAPE OF A DIFFUSER LOBE FOR A FOURIER DECOMPOSITION OF 
THE SHAPE INTO ITS FREQUENCY REPRESENTATION. 

FOREGO - COMPUTES THE COEFFICIENTS OF THE FOURIER SERIES 

REPRESENTATION OF R VS. ARC. ONE LOBE OF A DIFFUSER 
CROSS-SECTION IS THUS APPROXIMATED FROM PHI=-THETA TO 
PH 1= THETA. 

FORVAL - COMPUTES THE VALUE OF A FOURIER SERIES AT ARCS USING NTERM 
TERMS IN THE EXPANSION 

RWINPT - PERFORMS READING OR WRITING OF USER INPUT 

REPLIC - REPLICATES THE SINGLE LOBE REPRESENTATIVE INTO A MIXER 
CROSS-SECTIONAL VIEW 

OUT 30 - CREATES 3D OUTPUT FILE FOR 'MOVIE* PLOTTING 

LAMDAS - CALCULATES AND PLOTS UMBDAS AND LAMBDA PRIMES 


ROUTINES CALLED BY: 
FLOMIX PROGRAM 


RWINPT C AXITLT , XPCNT , STA , ISHAPE » ALF , THE , E Y , EZ , RE , FY , FZ , 
RF , IRWOPT ) 


THIS SUBROUTINE PERFORMS READING OR WRITING OF USER INPUT. 


INPUT AND OUTPUT ARGUMENTS: 

AXITLT, XPCNT, STA, ISHAPE, ALF, THE, EY.EZ, RE, FY,FZ,RF 
IRWOPT - READ /WRITE OPTION 
s 0, READ 
* 1, WRITE 


ROUTINES CALLEO BY: 

FORMIX - READS THE SPECIFIED DATA FILE, OR HELPS THE USER CREATE 
A NEW DATA FILE WHEN ONLY THE CONTOUR IS GIVEN. 


C************************************************************* 
SUBROUTINE 0UT30(XFF ,YFF ,ZFF ,PF,IANL,NSTA,NTOTL) 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CREATES 30 OUTPUT FILE FOR 'MOVIE 1 PLOTTING 

INPUT ARGUMENTS: 

XFF - X COORDINATES 

YFF - Y COORDINATES 

ZFF - Z COORDINATES 

PF - COMPLEX ARRAY FOR Y AND Z 

IANL - ARRAY IDENTIFYING STATIONS THAT WERE ANALYZED 

NSTA - NO. OF STATIONS 

NTOTL - NO. OF POINTS ABOUT A LOBE 

ROUTINES CALLED: 

REPLIC - REPLICATES THE SINGLE LOBE REPRESENTATIVE INTO A MIXER 
CROSS-SECTIONAL VIEW 


C ROUTINES CALLED BY: 
C FORMIX 



c ***** * * **♦#<»* 

SUBROUTINE AREA* X0*W,APRI > AF AN , AGAP , PEN , NLOBE . NST , IERR ) 

C ' 

C THIS SUBROUTINE CALCUUTES MIXER AREA ON A NON- VERTICAL PUNE 
C DEFINED BY THE X AXIS INTERCEPT , XO, AND ANGLE. W. 

C NOTE: POSITIVE U IS CCW FROM VERTICAL. 

C 


C INPUT ARGUMENTS: 

C XO - X AXIS INTERCEPT 
CM ANGLE 

C NLOBE - NO. OF LOBES 

C NST - NO. OF LOBE STATIONS 

C 

C OUTPUT ARGUMENTS: 

C APRI - PRIMARY AREA 

C AFAN - FAN VALLEY AREA 

C AGAP - GAP AREA 

C PEN - AMIX / ATOT 

C IERR - ERROR CODE (O-NO ERROR, 1 -ERROR) 

C 

C ROUTINES CALLED: 

C SPLIT - DETERMINES THE INTERSECTION OF A PUNE DEFINED BY 
C X-INTERCEPT, XO, AND ANGLE, U. , AND A SPLINE FIT CURVE 

C DEFINEO FROM VALUE OF INO. 

C 

C ROUTINES CALLED BY: 

C LOBGEN - CALCUUTES MIXER LOBE SHAPES AND CROSS-SECTIONS 
C 

ftWWWWwww ItwwW WWWwWWwWwwwW WWWWwWwWWWWWwW WWWwWWWwWwWW VfwWWWlVW fRWwW 


SUBROUTINE SPLIT* XO ,W, INO , XI ,R ,NST) 

c 

C THIS SUBROUTINE DETERMINES THE INTERSECTION OF A PUNE 
C DEFINE BY X INTERCIrT, XO, AND ANGLE, H. , AND A SPLINE 
C FIT CURVE DEFINED FROM VALUE OF INO. 

C IND IS USED IN FUNCTION RADIUS* X , IND ,R ) AS AN INDICATOR 
C OF WHICH SPLINE FIT IS TO BE READ. 

C 

C INPUT ARGUMENTS: 

C XO - X INTERCEPT 

C W ANGLE 

C IND - SPLINE FIT INDICATOR 

C NST - NO. OF LOBE STATIONS 

C 

C OUTPUT ARGUMENTS: 

C XI - X INTERSECTION 

C RI INTERSECTION RAOIUS 

C 

C ROUTINES CALLED: 

C RADIUS - EVALUATES SPLINE FIT CURVES AT SPECIFIED LOCATIONS 
C 

C ROUTINES CALLED BY: 

C AREA - CALCULATES MIXER AREA ON A NON-VERTICAL PLANE DEFINED 
C BY THE X AXIS INTERCEPT, XO, AND ANGLE, W. 

C************************************************************* 

REAL FUNCTION ZPRI CX.Y) 

c 

C THIS FUNCTION EVALUATES MIXER PRIMARY LOBE WIDTH AT AXIAL STATION X 
C AND HEIGHT Y. 

C USED AS F ( X ) FOR SIMPSON* S RULE INTEGRATION. 

C 

C INPUT ARGUMENTS: 

C X AXIAL STATION 

C Y HEIGHT 

C 

C OUTPUT ARGUMENTS: 

C ZPRI - MIXER PRIMARY LOBE WIDTH 
C 

C ROUTINES CALLED: 

C GEOM - CALCUUTES VERTICAL SECTION GEOMETRY 
C 

C ROUTINES CALLED BY: 

C AREA - CALCULATES MIXER AREA ON A NON-VERTICAL PUNE DEFINED BY 
C THE X AXIS INTERCEPT, XO, AND ANGLE, H. 

C 
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C ******************************************* ****************** 

SU BROUTINE LOBGENI STA.RXD , RE L .RFV.ROO , XPCNT ,NST , ITRANS » ALF »THET , 
BYY.BZZ.RBB.CYY.CZZ.RCC) 

C 

C RAOIAL HALL MIXER OECK - OBTAINED FROM D. HILEY <2-23-B3> 

C 

C THIS ROUTINE CALCULATES MIXER LOBE SHAPES AND CROSS-SECTIONS. 

C 

C INPUT ARGUMENTS: 

C STA - STATION X VALUES 

C RID - PLUS RADIUS 

C REL - ENGINE LOBE RADIUS 

C RFV - FAN VALLEY RADIUS 

C ROO OUTER CASE RADIUS 

C XPCNT - LOBE ANGLE PERCENTAGE 

C NST - NO. OF STATIONS 

C 

C OUTPUT ARGUMENTS: 

C ITRANS - MIXER SHAPE AT EACH AXIAL STATION 

C ALF - ANGLE TO MIXER HALL 

C THET - ANGLE OF HALF A LOBE 

C BYY - Y COORDINATE OF CENTER POINT OF LOBE 

C BZZ - Z COORDINATE OF CENTERPOINT OF LOBE 

C RBB - RAOIUS OF THE LOBE 

C CYY - Y COORDINATE OF CENTERPOINT OF FAN VALLEY 

C CZZ - Z COORDINATE OF CENTERPOINT OF FAN VALLEY 

C RCC - RADIUS OF FAN VALLEY 

C 

C ROUTINES CALLED: 

C BMFIT - SPLINE FITS MIXER LINES 
C GEOM - CALCULATES VERTICAL SECTION GEOMETRY 

C AREA - CALCULATES MIXER AREA ON A NON-VERTICAL PLANE 

C 

C ROUTINES CALLED BY: 

C FORMIX 


SUBROUTINE GEOM ( X ,STA, RID , REL, RFV, ROD .XPCNT, ISTA, ITRANS. IERR I 

c 

C THIS ROUTINES CALCULATES VERTICAL SECTION GEOMETRY. 

C 

C INPUT ARGUMENTS: 

C X - AXIAL STATION TO BE EXAMINED 

C STA - STATION X VALUES 

C RID - PLUG RADIUS 

C REL - ENGINE LOBE RADIUS 

C RFV - FAN VALLEY RAOIUS 

C ROD - OUTER CASE RADIUS 

C XPCNT - LOBE ANGLE PERCENTAGE 

C ISTA - NO. OF STATIONS 

C 

C OUTPUT ARGUMENTS: 

C ITRANS - MIXER SHAPE AT STATION X 
C IERR - ERROR INDICATOR 

C 

C ROUTINES CALLED: 

C BEVALE - EVALUATES BMFIT SPLINES 

C ARC - CALCULATES VALUES FOR LOBE RAOIUS. FAN VALLEY RAOIUS. AND 

C THE HEIGHT OF THE FAN VALLEY. 

C RAOIUS - EVALUATES SPLINE FIT CURVES AT SPECIFIED LOCATION 
C 

C ROUTINES CALLED BY: 

C LOBGEN - CALCULATES MIXER LOBE SHAPES AND CROSS-SECTIONS. 

C 
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SUBROUTINE ARC (RS.RFAN) 

C 

C GIVEN THE SIDE VIEW ENGINE LOSE HEIGHT, RE LI, AND A POINT 
C ON THE RADIAL WALL, RS, WHERE THE ENGINE LOBE RADIUS, RBI, 

C AND THE FAN VALLET RADIUS, RCI, MUST BE TANGENT, THIS ROUTINE 
C CALCULATES THE VALUES FOR RBI AND RCI, AND THE HEIGHT OF THE 
C FAN VALLEY, RFAN. 

C 

C INPUT ARGUMENTS: 

C RS POINT ON RADIAL WALL 
C 

C OUTPUT ARGUMENTS: 

C RFAN - HEIGHT OF FAN VALLEY 
C 

C ROUTINES CALLED: 

C 

C ROUTINES CALLED BY: 

C GEOM - CALCULATES VERTICAL SECTION GEOMETRY 
C 


*|AL_FWCTIW_R«IUS(XIN,INO,STA,RID,REL,RFV,ROD,XPCNT,ISTA) 

C THIS FUNCTION EVALUATES SPLINE FIT CURVES AT XIN. 

C I HO IS USED AS INDICATOR DEFINING CURVE IDENTITY 
C 

C INPUT ARGUMENTS : 

C XIN - X VALUE AT WHICH SPLINE FIT CURVES ARE EVALUATED 
C IND - INDICATOR DEFINING CURVE IDENTITY 

C STA - STATION X VALUES 

C RIO - PLUG RADIUS 

C REL - ENGINE LOBE RADIUS 

C RFV - FAN VALLEY RADIUS 

C ROD - CUTER CASE RADIUS 

C XPCNT - LOBE ANGLE PERCENTAGE 

C ISTA - NO. OF STATIONS 

C 

C ROUTINES CALLED: 

C BEVALE - EVALUATES BMFIT SPLINES 
C 

C ROUTINES CALLED BY: 

C GEOM - CALCULATES VERTICAL SECTION GEOMETRY 
C 

DIMENSION STAC 1 ),RID< 1 ),REL< 1 ),RFV( 1 ),RCD( 1 ). XPCNT ( 1 1 
COMMON /SPLIN2/ ARID ( 50 ) , BRIDC50), CPIDC50), DRIDC50I, 

2 ARELC50 ) , BRELI50), CRELI50I, DRELC50 ) , 

3 ARFVC50), BRFVCSO), CRFVC50I, ORFVC50), 

<► ARODC50), BRODCSO), CROOC50J, 0R0DC50), 

5 APCNTC 50 ) , BPCNTC50), CPCNTC50 ) , OPCNTC50 ) 

COMMON /SLOPE/ RIDIP, RELIP, RFVIP, RCDIP 


IND = 1 , READ PLUG RADIUS 

= 2, REAO ENGINE LOBE RADIUS 
= 3, REAO FAN VALLEY RADIUS 
= 4, REAO OUTER CASE RAOIUS 
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C#******«*##*#«***#*#*#*#*#****#****###******«****** i,,HHHHHHHH * 

SUBROUTINE EQANGS ( EOA»EAR »FOA»FAP ,THETAO»ALPHA,ISIDOP»CURRAO*N,P) 
C 

C THIS SUBROUTINE COMPUTES THE SERIES OF <Z,Y> PAIRS KWICH DESCRIBE THE 
C SHAPE OF A DIFFUSER LOBE FOR A FOURIER DECOMPOSITION OF THE SHAPE 
C INTO ITS FREQUENCY REPRESENTATION. 

C 

C INPUT; 

C EOA = 

C EAR = 

C FOA = 

C FAR = 

C THETAO = 

C ALPHA = 

C 

C ISIDOP = 

C CUR"* AO = 

C N = 

C 
C 

OUTPUT: 

P(I) = 


Y-DISTANCE TO CENTER OF CIRCULAR ARC AT THE PEAK, 

RADIUS OF THE CIRCULAR ARC AT THE PEAK, 

Y-DISTANCE TO CENTER OF CIRCULAR ARC AT THE VALLEY, 
RADIUS OF THE CIRCULAR ARC AT THE VALLEY, 

ANGLE FROM THE PEAK TO THE VALLEY, 

ANGLE TO THE END OF THE UPPER ARC AND THE START OF THE 
LOWER ARC 

SIDE LOBE GEOMETRY OPTION 

SIDE LOBE RADIUS IF CURVED GEOMETRY CISIDOP=2) 

TOTAL NUMBER OF POINTS FROM THE PEAK TO THE VALLEY EQUAL 
SPACED IN ARC LENGTH. 


( Y »Z ) PAIRS DEFINING ONE LOBE OF THE MIXER FROM 


ROUTINES CALLED: 

ROUTINES CALLED BY: 
FORMIX 


C««**«*««tf*«««***«*«««*«««*^«^*M*««****«*«**«************** 
SUBROUTINE LAHDAS C IANL ) 

C — — 

C CALCULATES AND PLOTS LAMBDAS AND UMBO A PRIMES 
C 

C INPUT ARGUMENTS: 

C IANL - ARRAY IDENTIFYING STATIONS THAT WERE ANALYZED 

C 

C ROUTINES CALLED: 

C BMFIT - PERFORMS A SPLINE FIT OF THE USER LAMBDAS 

C BEVAL - PERFORMS AN EVALUATION BASED ON THE SPLINE FIT 

C 

C ROUTINES CALLED BY: 

C FORMIX 
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SUBROUTINE FOREGO C NPl , P , ARSIN , BRCOS , AZR , B NR ) 

C “ 

C THIS SUBROUTINE COMPUTES THE COEFFICIENTS OF THE FOURIER SERIES 
C REPRESENTATION OF R VS. ARC. ONE LOBE OF A DIFFUSER CROSS-SECTION 
C IS THUS APPROXIMATED FROM PHI=-THETA TO PHI=THETA. 

C 

C INPUT: 

C NPL 

C P(I) 

c 

C OUTPUT: 

C ARSIN 

C BRCOS 

C 

C AZR 

C BNR 

C DELPHI 

C 

C ROUTINES 

C SERIES 

C 

C ROUTINES 

C FORMIX 

C 


c 

SUBROUTINE FORVAL ( NTERM, ANGL , ACOEF ,BCOEF , AZ , FVAL>BNR ) 

c 

C THIS SUBROUTINE COMPUTES THE VALUE OF A FOURIER SERIES AT ARCS 
C USING NTERM TERMS IN THE EXPANSION. 

C 

C INPUT: 

C NTERM = NUMBER OF TERMS USED IN THE EXPANSION. 

C ANGL = ANGLE TO BE EVALUATED 

C ACOEF ( I ) = COEFFICIENT OF COSINE TERMS, 1=1 , . . . ,NMAX. 

C eCOEF(I) = COEFFICIENT OF SINE TERMS, 1=1 , . . . >NMAX. 

C AZ = CONSTANT TERM IN THE SERIES. 

C 

C OUTPUT: 

C FVAL = VALUE OF THE EXPANSION. 

C 

C ROUTINES CALLED: 

C 

C ROUTINES CALLED BY: 

C FORMIX 

C 


= EVEN NUMBER OF OATA POINTS. 
« (Z,Y) PAIRS, 1=1 NTOTL. 


= ( NPL/2-1 ) COEFFICIENTS FOR THE SINE TERMS OF THE EXPANSION 
= (NPL/2-1) COEFFICIENTS FOR THE COSINE TERMS OF THE 
EXPANSION 

= CONSTANT TERM IN THE SERIES. 

= COEFFICIENT BNZ OF BNR*COS< NTOTL/2 "DELPHI ) , WHERE 
= 6.2S31B*(K-1 )/NPL» K=1 ,2 , . . . ,NPL. 

CALLED: 

- DETERMINATION OF COEFFICIENTS IN SERIES 
CALLED BY: 


C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 


SUBROUTINE REPLIC ( PF .NTOTL ,NLOBE , IP, I3D ,XF ) 

S " ,StE t0M INTO 


(ZF»YF) PAIRS, 1=1,... , NTOTL. 

STATION INDICATOR 

SINGLE OR TOTAL STATION REPLICATION INDICATOR 
NO. OF PAIRS 
STATION LOCATIONS 

ROUTINES CALLED: 

PLTEK2 

ROUTINES CALLED BT: 

FORMIX 


INPUT: 

PF( I ) = 
IP = 
1 3D = 
NTOTL = 
XF = 
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SUBROUTINE SERIES C Y»N. AZ* A*B* AN. YY *Z*M ) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DETERMINATION OF COEFFICIENTS IN SERIES 


INPUT: 
YIK ) 


N 


INPUT VALUE OF FUNCTION AT AN ANGLE OF TK=2 .*PI*( K-l )/N 
FOR K=1 TO K=N 
PI = 3.1415927 

EVEN NUMBER OF INPUT FUNCTION VAIUES (MUST EXCEED 2) 


OUTPUT 

AZ 

A 

B 

AN 

YY 

Z 

u 


CONSTANT TERM IN SERIES 

(N/2-1) REAL OUTPUT VALUES OF COEFFICIENTS IN C0S *??*** 
(N/2-1) REAL OUTPUT VALUES OF COEFFICIENTS IN SINE SERIES 
COSINE (N/2*TK) TERM IN SERIES 
DUMMY STORAGE OF LENGTH N*2 (REAL) 

DUMMY STORAGE OF LENGTH 2*N (REAL) 

U = DUttMY STORAGE OF LENGTH N (REAL) 

L=N/2-l 

Y(KI=AZ* SUM (ACL)«COStL«TK)*B<LM»SINCL*TKn*AN»COSCM/2*TK) 

*mminMHHH**#***i*******^ ******* *****^******^*^^^* , **^* < *** , * , * ,MH ** , ****** J **** 


ROUTINES CALLED: 

RFAST - FAST FOURIER TRANSFORM OF REAL DATA 


ROUTINES 

FOREGO 


CALLED BY: 

- COMPUTES THE COEFFICIENTS OF 
REPRESENTATION OF R VS. ARC. 
SECTION IS THUS APPROXIMATED 


THE FOURIER SERIES 
ONE LOBE OF A DIFFUSER CROSS 
FROM PHI=-THETA TO PHI=THETA. 


C 

C 

c 

C 

c 

c 

C 

c 

C 

c 

c 

C 

C 

C 

c 

C 

c 

c 

c 

C 

c 

c 

c 

c 


SUBROUTINE SERINV( Y »N> AZ> A»B* AN*YY *Z*W) 


REVERSE OF ‘SERIES’ 


Y(K) = INPUT VALUE OF FUNCTION AT AN ANGLE OF TK=2 .*PI*(K-1 )/N 
FOR K=l TO K=N 

PI = 3.1415927 _ _ 

N s EVEN NUMBER OF INPUT FUNCTION VALUES (MUST EXCEED 2) 


OUTPUT: 

AZ = CONSTANT TERM IN SERIES 

A = (N/2-1) REAL OUTPUT VALUES OF COEFFICIENTS IN COS SERIES 

B = (N/2-1) REAL OUTPUT VALUES OF COEFFICIENTS IN SINE SERIES 

AN = COSINE (N/2*TK) TERM IN SERIES 
YY = DUMMY STORAGE OF LENGTH N*2 (REAL) 

Z = DUMMY STORAGE OF LENGTH 2*N (REAL) 

H = DUMMY STORAGE OF LENGTH N (REAL 

ROUTINES CALLED: 

R FASTI - REVERSE OF ’RFAST* 

ROUTINES CALLED BY: 
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SUBROUTINE PFASTCX, Y,N,Z,W,S) 

C 1 

C FAST FOURIER TRANSFORM OF REAL DATA 
C 

C INPUT: 

C X = N REAL INPUT VALUES 

C N = EVEN NUMBER OF INPUT VALUES (MUST EXCEED 2) 

C S = SIGN CONTROLLING DIRECTION OF TRANSFORM 

C 

C OUTPUT: 

C Y = N/2*l COMPLEX OUTPUT VALUES 

C Z = DUMMY STORAGE OF LENGTH 2N (REAL) 

C W = DUMMY STORAGE OF LENGTH N (REAL) 

C 

C THIS PRODUCES 'OUTPUT Y' FROM 'INPUT X', WHERE 

C ******»**«****Ht*«********«**tt«*tt*«*««««**«*#««»*«*«iHtff«****«tf«**** 

C K=N 

C Y( J)=SUM X(K)*EXP(SIGN< 1. ,S)*I«2*PI»< J-l)*(K-l)/N) 

C K = i 

C *###*«*##«IH»*«#«#**#*»Ht#**##*******#*#»**#*###«#«#*#*##*# J HHH**#«*^ 

C WITH I=SQRT( -1 ) AND PI=3. 14159 

C NOTE THAT Y< N-Jf2 )=CONJ( Y( J )) FOR J=1 TO J=N/2U 
C THUS ONLY YU) TO Y(N/2*1> ARE CALCULATED 

C COMPLEX NUMBERS ARE HANDLED IN FORTRAN 4 CONVENTION, NAMELY THE 
C REAL AND IMAGINARY PARTS ARE STORED IN ALTERNATE CELLS, STARTING 
C WITH THE REAL PART OF Y(l) IN THE FIRST LOCATION, ETC. 

C 

C ROUTINES CALLED: 

C FAST - FAST FOURIER TRANSFORM OF COMPLEX DATA 
C 

C ROUTINES CALLED BY: 

C SERIES - DETERMINATION OF COEFFICIENTS IN A SERIES 
C 


C 

SUBROUTINE FAST (X,Y,N,2,W,S ) 

c 

C FAST FOURIER TRANSFORM OF COMPLEX DATA 
C 

C INPUT: 

C X = N INPUT VALUES (COMPLEX) 

C N = NUMBER OF VALUES 

C 

C OUTPUT: 

C Y * N OUTPUT VALUES (COMPLEX) 

C Z = DUMMY STORAGE OF LENGTH 2N (COMPLEX) 

C W = DUMMY STORAGE OF LENGTH N (COMPLEX) 

C S = SIGN CONTROLLING DIRECTION OF TRANSFORM 

C 

C THIS PR00UCE5 'OUTPUT Y' FROM 'INPUT X', WHERE 

C «****#***tftf**»tf**tf«**«ft**tf*****«*tf*«tf*tt«*«*«tf*tf«IMHHfffft***ft*««***** 

C K=N 

C Y( J )-SUM X(K)*EXP(SIGN( 1. ,S)*I*2*PI*< J-l )*< K-l )/N) 

C K=I 

C tf**«*«***tt**tf«*«tf**«***«*tf#«ffKtf«**««<Hf*«*ff»tf*M***ff ***#***»#****««» 

C WITH I=SQRT(-1), S=*l. OR S=-l., AND PI=3. 14159 

C 

C ROUTINES CALLED: 

C 

C ROUTINES CALLED BY: 

C RFAST - FAST FOURIER TRANSFORM OF REAL DATA 
C RFASTI - REVERSE OF 'RFAST' 

C 
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SUBROUTINE EBC 


CALCULATES THE POSITION OF THE DOWNSTREAM BOUNDARY AND THE DOWNSTREAM 
MACH NUMBERS FOR EACH STREAM. 

INPUT AERO DATA: 

PT1, PT2 > PINF, TT1, TT2, 6AM1, SAM2, Rl, R2 
OUTPUT ARGUMENTS: 

EM1S - PRIMARY STREAM MACH NO. FOR THE SPLITTER 
EM2S - SECONDARY STREAM MACH NO. FOR THE SPLITTER 

ROUTINES CALLED: 

HAF - SEARCH ROUTINE USED BY EBC 

ROUTINES CALLED BY: 

FLOMIX PROGRAM 


SUBROUTINE HAF(MAX.MIN.X,Y,IC.IER » 

c 

C THIS ROUTINE IS A SEARCH ROUTINE FOR EBC. 

C 

C INPUT ARGUMENTS: 

C MIN - MINIMUM VALUE OF OUTPUT PARAMETER FOR SEARCH 

C MAX - MAXIMUM VALUE OF OUTPUT PARAMETER FOR SEARCH 

C X CURRENT VALUE OF OUTPUT PARAMETER DURING SEARCH 

c Y INPUT PARAMETER TO BE MATCHED 

C IC - l» INPUT PARAMETER IS SPLITTER SECONDARY FLOW AREA 
C - 2, INPUT PARAMETER IS COMPOUND FLOW FUNCTION. BETA 

C - 3. INPUT PARAMETER IS SECONDARY FLOW RATE 

C 

C OUTPUT ARGUMENTS: 

C X - IC=1. OUTPUT PARAMETER IS SPLITTER STATIC PRESSURE 
C - IC-2. OUTPUT PARAMETER IS NOZZLE EXIT STATIC PRESSURE 

C - IC=3. OUTPUT PARAMETER IS DOWNSTREAM BOW® ARY STATIC PRESSURE 

C IER - ERROR INDICATOR: 0 IS O.K. ► 1 IS NO CONVERGENCE IN 20 TRIES. 

C 

C ROUTINES CAL LEO: 

C 

C ROUTINES CALLED BY: 

C EBC - CALCULATES THE POSITION OF THE DOWNSTREAM BOUNDARY AM) THE 
C DOWNSTREAM MACH NUMBERS FOR EACH STREAM 


C*««*«^•»**««*««***»********* ,, *** ,H, *******^^^^^^^^^^^^^^^^ 

SUBROUTINE GRID 

C THIS SUBROUTINE GENERATES A CARTESIAN MtSH FOR THE DENSEST GRID. 
C 

C ROUTINES CALLEO: 

C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
C 
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SUBROUTINE XPSUBS 
C 

C THIS SUBROUTINE 6ENERATES COARSE GRIO SUBSETS OF THE DENSEST 
C AS WELL AS CALCULATING COARSE GRID VALUES OF VARIABLES 
C SUCH AS IMAX, JMAX, IH, AND ITE. 

C 

c ROUTINES CALLED: 

C 

C ROUTINES CALLED BT: 

C FLOMIX PROGRAM 
C 


GRIDi 


SUBROUTINE BCSET3 1 XNTR .RNTR .IMXINT 1 
C ' 

C THIS SUBROUTINE GENERATES THE ARRAYS WHICH ARE USED TO ESTABLISH THE 
C EXACT SURFACE BOUNDARY COtCITIONS. 

C 


C OUTPUT ARGUMENTS: 

C XNTR - X INTERPOLATED SURFACE COORDINATES FOR CPPLOT 

C RNTR - R INTERPOLATED SURFACE COORDINATES FOR CPPLOT 

C IMXINT - NO. OF VALUES IN X AND R FOR EACH SURFACE 

C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
C 

C *•»•*»•»•»****«*«***«»«»•»*»«»«*.»»••••»•»««« 

C CONTOUR BOUNDARIES 

C L = I - CENTERBOOY 

C = 2 - INTERNAL SPLITTER 

C = 3 - EXTERNAL SPLITTER 

C = 4 - NOZZLE 

C*«**«*#tf *«»**«*«********«*»**«*«*«**««»«#«•** 

c 


SUBROUTINE PGUESS 

c 

C THIS SUBROUTINE GENERATES AN INITIAL GUESS FOR THE POTENTIAL, PHI. 
C THE GUESS IS OBTAINED EITHER FROM A SAVED SOLUTION, THE 
C SOLUTION FROM THE PREVIOUS GRID, OR A GIVEN FUNCTIONAL FORM. 

C 

C ROUTINES CALLED: 

C 

C ROUTINES CALLED BT: 

C FLOMIX PROGRAM 
C 


SUBROUTINE PLTGRDC XNTR , RNTR , IMXINT ) 

THIS ROUTINE HANDLES GRID PLOTTING 
INPUT ARGUMENTS: 

XNTR - X INTERPOLATED SURFACE COORDINATES FROM BCSETS 

RNTR - R INTERPOLATED SURFACE COORDINATES FROM BCSETS 

IMXINT - NO. OF VALUES IN X AND R FOR EACH SURFACE 

ROUTINES CALLED: PLTEK2 

ROUTINES CALLED BY: FLOMIX PROGRAM 
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SUBROUTINE INITL 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS SUBROUTINE SETS UP INITIALIZATION FOR EXTRA POUTED RELAXATION. 
ROUTINES CALLED: 


ROUTINES CALLED BY: 

FLC^IX PROGRAM 

1 IF AN EXTRAPOLATION DOES NOT OCCUR WITHIN INSTEP SWEEPS, AN 
EXTRAPOLATION WILL BE FORCED. AT THIS POINT A NEW TOLJMP 
COEFFICIENT WILL ALSO BE COMPUTED. 

2. IXTRAP - A SWEEP COUNTER FOR COMPARISON WITH INSTEP 

3. AT LEAST INOJMP SWEEPS WILL OCCU? BETWEEN EXTRAPOLATIONS. 

4 ! TOLJMP - A COEFFICIENT WHICH IS MULTIPLIED BT THE STANDARO 

DEVIATION OF THE DISTRIBUTION OF CORRECTIONS OVER THE MATRIX FOR 
A GIVEN SWEEP. 

THE ABS. VALUE OF THE AVERAGE CORRECTION IS COMPARED WITH TOLJMP 
TIMES THE STANDARD DEVIATION. 

NSTEPS - THE NUMBER OF NOOES IN THE GRID. 

SSTW™ OF THE RATIOS OF THE CORRECTIONS WITHIN EACH 
SWEEP . 

RLAMSS - THE SUM OF THE SQUARES OF THE RATIOS OF THE CORRECTIONS 

WITHIN ANT GIVEN SWEEP. 

"LAMS - -T4E S*r. U5 THE CORRECTS OVER ANY GIVEN SWEEP. 

ISTEP - A COUNTER WHICH IS COMPARED WITH INDJMP AND RESET TO I 
WHEN AN EXTRAPOLATION TAKES PLACE. 

WSWAVE - THE NUMBER OF SWEEPS USED IN CALCULATING THE AVERAGE 
RESIDUAL FOR AN EXTRAPOLATION. 

ISWAVE - A COUNTER FOR COMPARISON WITH NSWAVE. 


5 . 

6 . 

7 . 

A. 

9 . 

U. 

11 . 


SUBROUTINE SLOP 
C 

C THIS SUBROUTINE PERFORMS A SUCCESSIVE LINE OVER-RELAXATION. FOR 
C EACH COLUMN IN THE CARTESIAN MESH, A TRIDIAGONAL MATRIX 
C IS SET UP AND SOLVED. THE CORRECTIONS ARE APPLIED WITH 
C AN OVER-RELAXATION FACTOR FOR SUBSONIC FLOW, AND AN 
C UNDER-RELAXATION FACTOR FOR SUPERSONIC FLOW. WHEN 
C THE PROPER CRITERIA ARE MET, AN EXTRAPOLATION IS CALLED. 

C 

C ROUTINES CALLED: 

C MATRIX - SETS UP THE MATRIX OF COEFFICIENTS RESULTING 
C FROM THE PARTIAL DIFFERENTIAL EQUATION AM) THE BOUNDARY 

C CONDITIONS. THE R*PHIR RULE FOR THE CENTERBOOY AND THE 

C OUTER BOUNDARY IS INCLUDED IN THE MATRIX TO SIMPLIFY THE 

C LOGIC FOR VARIABLE COLUMN LENGTHS. AT SOME FUTURE DATE 

C IT MAY BE CHANGED FOR THE SAKE OF EFFICIENCY. 

C TRIDI - SOLVES A SET OF N TRIDIAGONAL EQUATIONS AS OBTAINED IN 
C THE MAIN PROGRAM FOR A RADIAL COLUMN. 

C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
C 


C**»* *****»»«»*«**»«<•»»*«»«»****«***«»» »*»******»»«»**«** 

SUBROUTINE FARFLD 
C 

C THIS SUBROUTINE GENERATES THE FAR FIELD BOUNDARY CONDITIONS. IT IS 
C INCLUDED IN THE RELAXATION SWEEPING LOOP TO ALLOW FOR PERIOOIC 
C UPDATING BASED ON THE CIRCULATION. 

C 

C ROUTINES CALLED: 

C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
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SUBROUTINE MATRIX 
C 

C THIS SUBROUTINE SETS UP THE MATRIX OF COEFFICIENTS RESULTING 
C FROM THE PARTIAL DIFFERENTIAL EQUATION AND THE BOUNDARY 
C CONDITIONS. THE R*PHIR RULE FOR THE CENTER BODY AND THE 
C OUTER BOUNDARY IS INCLUDED IN THE MATRIX TO SIMPLIFY THE 
C LOGIC FOR VARIABLE COLUMN LENGTHS. AT SOME FUTURE DATE 
C IT MAY BE CHANGED FOR THE SAKE OF EFFICIENCY. 

C 

C ROUTINES CALLED: 

C BCABCM - ESTABLISHES BOUNDARY CONDITION VALUES FOR MATRIX 
C COEFFICIENTS A, B» C, AND H 

C POEABC - CALCULATES THE POE MATRIX COEFFICIENTS A,B, AC, AND THE 
C RIGHT HAND SIDE W 

C 

C ROUTINES CALLED BY: 

C SLOR - PERFORMS A SUCCESSIVE LINE OVER-RELAXATION. FOR 
C EACH COLUMN IN THE CARTESIAN MESH, A TRIDIAGONAL MATRIX 

C IS SET UP AND SOLVED. THE CORRECTIONS ARE APPLIED WITH 

C AN OVER-RELAXATION FACTOR FOR SUBSONIC FLOW, AND AN 

C UNDER -RELAXATION FACTOR FOR SUPERSONIC FLOW. WHEN 

C THE PROPER CRITERIA ARE MET, AN EXTRAPOLATION IS CALLED. 

SUBROUTINE TRIDI C A, B,C, VECTOR ,H,N ) 

c 

C THIS SUBROUTINE SOLVES A SET OF N TE. DIAGONAL EQUATIONS AS 
C OBTAINED IN THE MAIN PROGRAM FOR A RADIAL COLUMN. THE RESULTS 
C ARE LEFT IN THE ARRAY »VECTOR* AND A,B,C AND W ARE DESTROYED. 

C 

C INPUT: 

C A - WORK VECTOR FOR TRIO I AGONAL EQUATIONS 

C B - WORK VECTOR FOR TRIDIAGONAL EQUATIONS 

C C - WORK VECTOR FOR TRIDIAGONAL EQUATIONS 

C W - WORK VECTOR FOR TRIDIAGONAL EQUATIONS 

C N - NO. OF TRIDIAGONAL EQUATIONS 

C 

C OUTPUT: 

C VECTOR - RESULTING ARRAY 
C 

C ROUTINES CALLED: 

C 

C ROUTINES CALLED BY: 

C SLOR - PERFORMS A SUCCESSIVE LINE OVER-RELAXATION. FOR 
C EACH COLWN IN THE CARTESIAN MESH, A TRIDIAGONAL MATRIX 

C IS SET UP AND SOLVED. THE CORRECTIONS APE APPLIED WITH 

C AN OVER-RELAXATION FACTOR FOR SUBSONIC FLOW, AND AN 

C UNDER-RELAXATION FACTOR FOR SUPERSONIC FLOW. WHEN 

C THE PROPER CRITERIA ARE MET, AN EXTRAPOLATION IS CALLED. 

************ a#*****#****^**##****##,!^*#**##**#******* 

SUBROUTINE RESCL1 


THIS SUBROUTINE CALCULATES RESIDUES FOR EACH SWEEP THROUGH RELAXATION 
LOOP. 

ROUTINES CALLED: 

MATRIX - SETS UP THE MATRIX OF COEFFICIENTS RESULTING 

FROM THE PARTIAL DIFFERENTIAL EQUATION AND THE BOUNDARY 
CONDITIONS . THE R*PHIR RULE FOR THE CENTERBOOY AND THE 
OUTER BOUNDARY IS INCLUDED IN THE MATRIX TO SIMPLIFY THE 
LOGIC FOR VARIABLE COLUMN LENGTHS. AT SOME FUTURE DATE 
IT MAY BE CHANGED FOR THE SAKE OF EFFICIENCY. 

ROUTINES CALLED BY: 

FLOMIX PROGRAM 
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SUBROUTINE PD E ABC! ISTR I 




THIS SUBROUTINE CALCULATES THE PARTIAL DIFFERENTIAL EQUATION 
MATRIX COEFFICIENTS* A*B» AM) C> AND THE 
RIGHT HAND SIOE. H. 


ROUTINES CALLED: 


MATRIX - SETS UP THE MATRIX OF COEFFICIENTS RESULTING 

FROM THE PARTIAL DIFFERENTIAL EQUATION AND THE BOUICARY 
CONDITIONS. THE R*PHIR RULE FOR THE CENTERBGDV AND THE 
OUTER BOUIOARY IS INCLUDED IN THE MATRIX TO SIMPLIFY THE 
LOGIC FOR VARIABLE COLUM LENGTHS. AT SOME FUTURE DATE 
IT MAY BE CHANGED FOR THE SAKE OF EFFICIENCY. 


»****»«»«»**«»**«****«»«*«***»•*•»»**•*«*« »wmhhhhhhmhmmhhhm» 

SUBROUTINE BCABCNt I5TR , IBNO ) 

THIS SUBROUTINE ESTABLISHES BOUM)ARY CONDITION VALUES FOR MATRIX 
COEFFICIENTS A* B. C. AM) N 

INPUT ARGUMENTS: 

JSTR - STREAM NUMBER t 1 -PRIMARY, 2-SECONDARY) 

IBM) - BOOM) ARY NUMBER ( X-CENT. ,2-INNEfi SPL. .3-OUTER SPL. , A- NOZZLE ) 
ROUTINES CALLEO: 

DELTAF - CALCULATES DELTA FUNCTIONS 
ROUTINES CALLEO BY: 

MATRIX - SETS UP THE MATRIX OF COEFFICIENTS RESULTING 

FROM THE PARTIAL DIFFERENTIAL EQUATION AM) THE BOUNDARY 
COYCITIONS. THE R*PHIR RULE FOR THE CENTERBOOY AND THE 
OUTER BOUNDARY IS INCLUDED IN THE MATRIX TO SIMPLIFY THE 
L06IC FOR VARIABLE COLUFM LENGTHS. AT SOME FUTLDE DATE 
IT MAY BE CHANGED FOR THE SAKE OF EFFICIENCY. 

l»*«ll«*«*«*««»*^•*««»«*«**•»^»»«^<«*■*•*•* ,, »** ,, *** ,,,, **** ,H, ****' MM, 

SUBROUTINE PLUME 


THIS SUBROUTINE IS USED FOR CONFIGURING THE PLIR1E. 

ROUTINES CALLEO: 

ROUTINES CALLEO BY: 

MATRIX - SETS UP THE MATRIX OF COEFFICIENTS RESULTING 

FROM THE PARTIAL DIFFERENTIAL EQUATION AND THE BOLAOARY 
CONDITIONS. THE R*PHIR RULE FOR THE CENTERBOOY AND THE 
OUTER BOUM>ARY IS INCLUDED IN THE MATRIX TO SIMPLIFY THE 
LOGIC FOR VARIABLE COLUtM LENGTHS. AT SOME FUTURE DATE 
IT MAY BE CHANGED FOR THE SAKE OF EFFICIENCY. 


SUBROUTINE UPDATE 
C ■ — 

C THIS SUBROUTINE CALCULATES THE CORRECTIONS FOR THE NEXT SWEEP. 

C 

C ROUTINES CALLED: 

C EXTRAP - CALCULATES THE EXTRAPOLATED RELAXATION FOR THE NEXT SWEEP 
C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
C 


SUBROUTINE EXTRAP 

C THIS SUBROUTINE CALCULATES THE EXTRAPOLATED RELAXATION FOR THE NEXT 
C SUEEP. 

C 

C ROUTINES CALLED: 

C 

C # °UPDATE - CALCULATES THE CORRECTIONS LOR THE NEXT SUEEP 


C************************************************************* 
FUNCTION DELTAF ( I , J ) 

C "" r 

C THIS FUNCTIONS IS USED BY BCABCW TO CALCULATE DELTA FUNCTIONS. 
C 

C ROUTINES CALLED BY: 

BCABCW - ESTABLISHES BOtMDARY CONDITION VALUES FOR MATRIX 
COEFFICIENTS A, B, C, AND M 


C 

C 

c 

c* 


SUBROUTINE OUT PUT ( XNTR »RNTR , IMXINT 1 

c 

C THIS SUBROUTINE CALLS THE VARIOUS OUTPUT SUBROUTINES THAT APPLY TO 
C EACH GRIO. 

C 

C INPUT ARGUMENTS: 

C XNTR - X INTERPOLATED SURFACE COORDINATES FROM BCSETS 

C RNTR - R INTERPOLATED SURFACE COORDINATES FROM BCSETS 

C IMXINT - NO. OF VALUES IN X AND R FOR EACH SURFACE 

C 

C ROUTINES CALLED: 

C MASFLO - PERFORMS THE MASS FLOW INTEGRATIONS 

C CPPLOT - PLOTS CP CALCULATIONS 

C FLDPLT - PLOTS CPS IN THE FLOW AREA 

C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 
C 


C********************************************** tttf **«#»«**** 

SUBROUTINE PHIPLT ( PHI,X,R.IMAX,JMAX,NI,NJ.NK.K»IC.D,H.UP,SHR.ETA) 
C 

C THIS ROUTINE HANDLES PLOTTING OF PHI VALUES. 

C 

C ROUTINES CALLED: PLTEK2 

C 

C ROUTINES CALLED BY: 

C 


C**************** ****************************************** 
SUBROUTINE CPPLOT C XNTR , RNTR , IMXINT ) 

C 1 

C THIS ROUTINE PLOTS CP CALCULATIONS 
C 

C INPUT ARGUMENTS: 

C XNTR - X INTERPOLATED SURFACE COORDINATES FROM BCSETS 

C RNTR - R INTERPOLATED SURFACE COORDINATES FROM BCSETS 

C IMXINT - NO. OF VALUES IN X AND R FOR EACH SURFACE 

C 

C ROUTINES CALLED: PLTEK2 » MVCHAR 

C 

C ROUTINES CALLED BY: 

C OUTPUT - CALLS THE VARIOUS OUTPUT SUBROUTINES THAT APPLY 

C TO EACH GRID. THESE INCLUDE THE CP PLOTS, THE MASS 

C FLOW INTEGRATIONS, AND THE DRAG CALCULATION. 

C 


C********************************************************** 

SUBROUTINE FLDPLT 

c 

C PLOTS CPS IN THE FLOW AREA 
C 

C ROUTINES CALLED BY: 

C OUTPUT - CALLS THE VARIOUS OUTPUT SUBROUTINES THAT APPLY TO EACH 
C GRID 
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C»*«**»#****«*W***«*«*tt4MI*«««*tf*tf*«^tfft**«*«*******««*IHHHHHH»** 

SUBROUTINE RECOMB ( THE .CPS, IKXT.NT, NEXT ) 

C ' 

C THIS ROUTINE COMBINES EACH HARMONIC SOLUTION INTO ONE SOLUTION. 
C 

C INPUT ARGUMENTS: 

C THE - THETA ANGLES 

C CPS - CPS FOR EACH CONTOUR 

C I MX I NT - NO. OF VALUES IN EACH CONTOUR 

C NEXT - NO. OF EXTERNAL SPLITTER STATIONS 

C 

C ROUTINES CALLED BY: 

C FLOMIX PROGRAM 

r 


C*#*******#*************«*#*»**#**#*«************************ 

SUBROUTINE PJMPLT ( PJMPTT ,RMXPTT ,RAVPTT , IMXRSP , JMXRSP ,NSTOT »NHAPM ) 
C 

C THIS ROUTINE HANDLES PLOTTING OF THE JUMP IN PHI VALUES AND THE 
C MAXIMUM AND AVERAGE RESIDUALS 
C 

C INPUT ARGUMENTS: 

C PJMPTT - JUMP IN PHI VALUES 

C RMXPTT - MAXIMUM RESIDUALS 

C RAVPTT - AVERAGE RESIOUALS 

C IMXRSP - MAXIMUM RESIOUAL POINTER FOR EACH SWEEP FOR I 

C JMXRSP - MAXIMUM RESIDUAL FOINTER FOR EACH SWEEP FOR J 

C NS TOT * MAXIMUM NO. OF POINTS TO PLOT 

C I !HARM - NO. OF HARMONICS ( MOOES ) 

C 

C ROUTINES CALLED: PLTEK2 

C 

C ROUTINES CALLED BY: FLOMIX PROGRAM 

C 
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APPENDIX 


Scarf Angle Analysis 

Mixing performance has been found to depend on the geometry character! sties of 
the lobe trailing edge plane. In particular, cutback or scarfing can be used 
to optimize mixer performance. Typically , the scarf angle (?) can vary + 15° 
from a radial cut. Analysis of such configurations is complicated by the 
irregular radial surface presented when viewed in a cylindrical frame of 
reference. An alternate representation, using a skewed system aligned to the 
lobe trailing edge avoids this problem. 



The governing equations can be transformed from the physical (x,r) to the 
computational plane (e,^) by means of the following transformation. 


£ = x - A r >i=l 

c 


(A. 1 ) 


where A =constant and c=c©s?= constant. The local velocities can 

then be related to these coordinates through chain rule differentiation. 


^ = ( 

3x 3x 


} + 


9vj 3_ ( 
ax 3rj 


) = 

d 1 


(A. 2a) 


or 


a ( ) = if 3 ( 

3-r 


) + a ( 

dr 9-r] 


) =-&.( ) 
4-JL 2 ( 

C. < 5^1 


(A. 2b) 

) 


p3* - t 3r=- A PViJl (A ' 3) 
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The term has been introduced to reflect the comprisable term from the 
s.d.t. expansion. Applying these flux definitions to a balance over a general 
element, in a cylindrical mesh, yield* 



Discriticizing reduces (A. 4) to 









(A. 5) 


where g x , g r are differenced along the mesh when expressed in terms of 
(A. 3). Similarly one can treat boundary intersecting elements in this skewed 
system. Recalling equation (22) 


where 



F-W 



N-b 



= O 



x ' k + 


is the "known" surface flux along the main surface Rm (x). The flux balance in 
the skewed system becomes 



- Rm AX 



- o 


where the local are as again represent only the exposed portion of the cell. 
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T t is desirable, in general, to have a skewed radial mesh over the entire 
calculation domain but boundary conditions are typically defined along radial 
upstream and downstream planes. The skewing transformation (A.l) must consider 
A=yVx,r). A generalization of the previous analysis is simplified if we 
introduce the flux balance in terms of the contravariant velocity components; 
those defined normal to the (E,*}) coordinates. Following the approach of 
Doria, equations (A. 2) applied to (x,r) yields four separate equation systems 
for the transformation metrics; eq. 


( ) -x 


l ~ <5* + v )x **) 

® ~^r X$ + *\r 



-.-Hr.- ,1 

?A-Vr 


If A= constant, these equations yield the following metrics relations 


I Xrj= s r f = o c 

Associated with this transformation are the Jacobian 


(A. 6) 


T =■ C. “ = C 


A. 7a) 


and the three invariant components of the metric tensor 


Q (i — Xj + — I 

a il « x £ ^y] + r 6 r ^ =° 

a 2Z - ^ - I (A. 7b) 

The conservation flux balance, in terms of the contravariant velocities (U,V) 
associated with (E,’"]) directions, is now given by the following components 


Ju A 5 * A € 4 ^ + A e jj^ f -Ag, 


TVA^j =■ ^ ^ A*] 


a "3^=- r f Ar i + A l[ 


(A. 8) 


The leading term in each expression reflects the contribution of the 
freestream velocity in the x-aligned perturbation potential formulation. One 
can compare (A. 8) and (A. 5) to check for consistency of the models. 


| 



-A 


+ c h 


A* 









( J* u Af ) ^ ^ * 


ay - a 3i 


w 


In the above analysis the freestream component is neglected. For a general 
element it will cancel out with the opposing face's contribution. Special 
consideration is needed for the surface intersecting case. 

The previous analysis can now be extended to include problems to wherein the 
mesh skewing varies. 


f - > - A (x,r i r 



r 

CCx.r) 


The metric invariants and Jacobian 

t „ = i-r U 

c 3x 

*!* - A 3J 

1 c JX 


(A.7a,b) can be reevaluated in terms of 

£ r = - A -jr SJ 

C z 3r 


vi r - 1 4 A 

c c lr 


69 



The exact relationship for t; , C, A can be obtained as follows. Consider the 
mesh generated from the lobe trailing edge to some downstream phase. Generate 
mesh skewing from fan valley contour. 



The equation of the lobe trailing edge is 

^ r ) =■ Xo(<V) + (r-r 0 ) A‘‘ 

Similarly the downstream mesh boundary is given by 


Xj Cr ) = X<j 


Linearly adjusting between these boundaries 



X, U) - X°(r 0 ) 
- y 0 CQ 


one obtains a general function expression for the x grid variation 

y(l) ' (>-X)Xo( r ) 

In this frame A , at any yj line, is given by 

All) = A (x r) — 1<--C, ■,) 
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